


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1991 


Advanced instrumentation concepts and their 
application to nuclear power plants. 


Soule, Ralph T. 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/25655 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


j ۹ 0 U DLEY research materials and institutional publications created by the NPS community. 
FW : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


۱۱ KNOX appointed — and published — scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http: //wwwenps.edu/library 









wid, Bats 
E ) A Pn ; d t r 4 ig ۷: 1 ور‎ Vei HN 
u.» 


1 
Cr 
0 


DE 


Li r3 P 
TOME T NE ^ uv: ZO 
A AAC 
o ME ۱ NETUS. May UNA LXI LE s 
3 Y i P 2 n PI rien, a Aa E ur. 
=t A [] Na, 1 DES u امه‎ FY A 
S 0 7n ۲ si FR" e ER rer ۷ 3: Au de 
“o 1, ايو‎ | o Pd E bl joi E Ae $ ZO FOE 
Drud Ak ANA AA A 
A 0 "I f P Ds sr 
v Pet PEE نف‎ x NS ۸: 
ae 1 کر ا‎ E ا‎ y 
Va ` d A yo * CO ^D E 4 € Nu oe RAS 
a LI ' Us i ۱ 4 P. a ۳ TL S C br En ^ 
1 MEE ST DOTNET m 
er Pe Aaa e A y 
0 















E 
^ 
scs 

= 

a‏ دا اپ اع 


b 


Lr 












4 


LA E i 


m "m 
LA "3 ab p.93. n * ^ Ls Yi 4 


€. 


0 
m] 







^» 
= ¢* 
aM 
RA 
e. a a 
x Vor 
LN JC 
PN 


i 
۳ r] 


Alan ره‎ 



































e E d E y 1 ag a" A 1 
0 A 4 د‎ 'Á qu SC "TE f} iy x A H ME -s E T A ۶ 
E "m ۶ & A. a ul Aa Py Ty E n ۱ P. T 
۱ 4 E] y sa” Je TY Py EE i ope ا‎ Fr Tu y " n 4 
r ^ 4 LR GE m LL Fa un Ai 5 تم‎ "on 
i x i à 4 Ains y” C A i y 1 à ۱ 
n] تایه سا‎ y ye ” F 1 A F a ار‎ t 
۱ k ba Gee ۰» یش ام‎ "Vy Hess 
۰ . 1 A Iud 0 LF A d 
, a g 3 E P 1 impe ید اج‎ 22 XU 1 OF 
m 0 و‎ n 4 , >. 2 Za E P EY | er 
y. L L " ac a 
۰ à f UE AE dur TO RT REP 
am o $8 ET N a 
E 4 E MAT] ^ : Ze a 
^ ۰" f ., ۱ تب‎ Au IE | À AS LES 
۱ P ý E i PO PT, P 
x یچ ی‎ ۸.۵ ge 1 P w ۹ a y», A y, NET el 4 Ae Poo 
LEE] > d Au D L] Ex kd x Bg PEL "EE wo IL EE RER betes AT ار‎ Lo NS 
By de at , g Y PA? LL ee ae © oe) TA cor eo n dab 
لا‎ , , E a d «9 ff de a ۹ $ 
۰ ie gs O ۲ او کي‎ ۱ 
t t %s ۴ LA 
۳ 8 P di 
m é = y m 2 5 >, 
LI [EI ft 
او‎ ee Eu 
7 A as P] 
1 NT m 1 : sr 
pos eier z 4d S j 3 B aes و ی‎ S A V 
n yis A 7 H t a e E O E m e 24.5 | e £ n ی‎ D Er 
A » J PD" ۳ ر چ‎ por am 17 
s ۳ "D = 60. . ` ر‎ ۰ | M d hy T 
4 ™ ., rn y E 
۲ ۳ EE D E TS E 
a g Pa) = é ۲ P P K a ون(‎ =: fi ie 
a ۴ PX « f Eid a e i ALLE. dí Very M à v 
M n c . 1 TAete oe. S. tay a tins core 
L Iu o EE s ۱ Tn pe da XA سر‎ epu 
(6 E A i ns E A =. Ln FAP ns ] Ca 2 1 


PIED "1 nr. me و‎ 9 ۳ ۰ 
2 T h ٤ > TEER TE Fer ی‎ Papen مت‎ p BER a 

































































. * 5 d ad ud Jis tod Saa ^ S Am he, 
4 EI ۰ Y y ۷4 t : ^. p i 7 3 ets >. TU hi lud 5 aid" 
"E یوب‎ AT EE. L^ ray "TEE e 
H E لوا‎ b L^ V ‘ MUR AI zii. E PAS AR S Áo ka 2] 
۰ ده‎ " % xp. nt ۳۳ E . la ye ji pr ade S 
E $i an: d nene Lar al sg LC ND aa Lr eod A er 
y D 4 ED (oy aed oom o A kn pae 
۰ D o" 0. oa hee ws تخر ما‎ ۲ 
۰ 1,3 بو‎ SIRO Quas ا ی ن ی‎ ih Me r E H m 
۰ ۲ T Lad A "ae aa ha A ve A ae, AS CIE 
^ g ۲ ۱ "ET E 1 wt A ۳ " " we. ve JA Pia a Ts "te 
5 ۰ ” «e mira EO u Mew y 
"n 4 "P me a . $ ۰ = d A. P 
1 TE یف‎ T cr ad e 5 EST ۰ m: یی‎ TRAE r4 PA 3 
2 £ P a Lb M pi f > o 69. P Fb P y ی‎ 
۳ ; B ag AL ' a y g E A rua - 1 we 3 مه‎ s Los ^ à VA e c 
> . e we 1 £g ید‎ : 9 x E A کر‎ w. LA Fe, $7 NL مش‎ a us 
0 t 0 E "IA NE de £ ۰ ۹ Tu nga VEF at or a a N E 
es ul PETI PEE M | ۳ N | Dd 2 ol 3 این هر(‎ 
mr T) T A ee 1 E RES Fab Mo UP PLE e m HM > ur 
enn sf SEE ee, et ix aene A eii ira 
..>a y اه‎ 1 . ol | viel. [NE "TIME LE ng + s m 
o. $ O ÓN dr FOU Y Wu 1 T مه‎ rp Pont p 
E v ve SL. a > bI «d ۲ ۳ هر‎ wi dado 
E. ^u ET em € y ی‎ ۳ "puc: Vac nee ES 
T le york r E srí ر‎ O o Vd ار‎ De pire ره‎ a = 
{ ۳ ۰ a 1 Irà ار‎ o A م‎ ft OP ۳ JU d 297 i MACC D dal s vert midi تیا‎ 
i d T. o E A AS و‎ DIST u PEE E proa e =) 
اب‎ oit LP PET a Ls رکچ‎ EEL i E porn a E: 
f AE 1 Louis s Y t DEN EP hf rn 4 rege t2 00, CAE 
f DES c وکن .اهم و و ویس‎ EA Ty 
REI E P d ah. -è A : ET 
E n a PS دی‎ EL T? A عمش‎ ash Mi a hol teh n de 
- 5 E [172 vs pt [a E d A 
ie vi " lal i n n oda CATA Pe 
0 Fe ۰ PP ^ اب‎ . y» Mo, 
Vit! LJ ۳۹ P AY " Aue E Ao. 
-a E! ۱۳ بط نله‎ e. Aa 2. ssi بو‎ E > PN = ^ vy hans وج‎ 
i i "T 3 rt AS ; رل‎ "d " Mi ar ات دون‎ £ ne ee e 2 LM poi. 
۰ 1 pv "a ar T^^. it مان‎ had و وا‎ IM TU ER 
AL a E UE pè o Fé ri»! pat EVA culos ^ END Ted مد‎ a My 
E 1%, 0 و‎ ٩ و‎ hcec دا‎ COR TEL TET UM ir cin de PE ue” TUN 
EE کم ر‎ S, A xo uf ¡9 19 e at ER E را ا‎ 
A 0 » . a Gey Ath MH nw wei gy Aine 07, سونو‎ te تا‎ gr C P RR 2 
à vy. wv feos a mr 3 $e 4v ion La dA PANA ern A eL en 
v a Jn DC ^ LUE SL nace a ey e 
5 bi . s "3 "E UE ad rn €. e itio L e 7 p AN y» 
> . A B^ A a P m a ui j m wer am Da bu. m: 
uL ,و‎ » e» E ۵ را‎ se = vy نان‎ ado a 0d at 
1 t5 » + h A ۲ . 4 -— v nv ARE a Pr EU 
p . 3 ES OS uch un or E A ^ 4 : 
g ۱ D 9 er wei: ee 
E ^ و‎ > AQ Ch LN DD e unt! ee A B: VC A poit SINE ول‎ 
1 . à reu. JA e ns BI tage 26 
E. 2 i p DUCATI RUOTE I nA 
Siete aed u 
TAS Ue : 
a 
سین‎ gan 
' b y 
۰ LI 
و‎ 
u 
T. r i = eb ean eka p 
3 A ` e a RT r ; r ا‎ ITO ERS ARE TL AE یه‎ ba al “wt 
ie Y تا‎ À T a da. 2 x دی‎ PEU EC pex dar ER EEE TEN Ca 
E < ¢ EN mae bee i Esperas í مه مه هچ جر ره‎ 
امد هم‎ 0 wm Kx 
da MEX Li MURS. rore 
I Rma du PE هو رن‎ VECES A 
LAS E ado m E 
AA ۱ 


het: 
| EI. A e LATI m zu EE RR AE 
a 1 et = 1 i = K Pae m PL ی ییو‎ ۳ et : 
x P , ^ P " ۳ s o qa. 4 ۱ 
m 1 ur زب‎ zx od p. 3 Y PAIS ATE 2 nee $ vs 
; e s H UD - Cot ee 
JO r ham Rn a Br o x. e & E i 


v0 


v 


. 
Lj 0 
P 
1 
pi 
D 
O 
P 
0 
0 
P 
۰ 
E 
" 1 
0 ' 
0 
1 
1 
۰ 
t 
0 
1 
* 
FA 
t A 
E 
1 
0 
0 
b 
0 
E 
٤ 0 
۰ 
۱ 0 
a 
۹ 
D 
. 
PES 
. 
۰ 
" ۰ 
E 
نت‎ É 
P 
Li 
P 
g 
E 
E M 
~*~ 
b 


wer 


kod 


Soe 


1 











ADVANCED INSTRUMENTATION CONCEPTS AND THEIR APPLICATION TO 
NUCLEAR POWER PLANTS 


by 
RALPH T. SOULE 


B.S. Electrical Engineering, U.S Naval Academy 
(1982) 
SUBMITTED TO THE DEPARTMENTS OF 
NUCLEAR ENGINEERING 
and 
ELECTRICAL ENGINEERING AND COMPUTER SCIENCE 


IN PARTIAL FULFILLMENT OF THE REQUIREMENTS 
FOR THE DEGREES OF 


MASTER OF SCIENCE IN NUCLEAR ENGINEERING 
and 
MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 


at the 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


May 1991 


© Ralph T. Soule, 1991. 
All rights reserved 


The author hereby grants to MIT permission to reproduce and 
to distribute copies of this thesis document in whole or in 


er TERN. 





ADVANCED INSTRUMENTATION CONCEPTS AND THEIR APPLICATION TO 
NUCLEAR POWER PLANTS 
by 
RALPH T. SOULE 


Submitted to the Departments of Nuclear Engineering and 
Electrical Engineering And Computer Science in partial 
fulfillment of the requirements for the degrees of Master of 
Science in Nuclear Engineering and Master of Science in 
Electrical Engineering 


ABSTRACT 


The characteristics of computer based instrumentation, or 
smart instruments, are investigated. Computer based 
instruments are distinguished by their ability to include a 
more complex model of the physical processes influencing the 
desired measurement than is possible with conventional 
instrument. Smart instrumentation is described with 
emphasis on illustrating its ability to improve data 
collection, storage, display, and evaluation. The 
application considerations of redundancy, consistency, 
noise, and filtering are also addressed. 


As an application example, a smart instrument for measuring 
Steam generator water level in a pressurized water reactor 
is designed. A model, accounting for the important 
processes affecting level measurement is developed. An 
error exists in the computation the fluid shear stresses, 
but the model calculations remain illustrative of those 
pertinent to smart instrument design. The model is stable 
for both steady state and transient conditions, but there 
are restrictions on the rate of the transient. Simulated 
level data is used to compare a simplified level instrument 
with the smart level instrument. The smart instrument is 
more accurate, but not by more than one percent. Methods 
that could be used by a smart instrument to recover from 
operation outside its model assumptions are discussed. 
Finally, recommendations are made for future work. 
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Introduction 


Chapter 1 


The objective of this research is to suggest a method for 
the application of microcomputer based instrumentation to 
nuclear power plants. The specific instrument illustrated 
is for a steam generator level measurement, but the basic 
Steps of physical process selection, variable 
identification, model creation, and integration of the 
process model with signal validation are applicable to any 
important plant parameter. A procedure to implement a smart 
instrument is suggested. A model was developed to include 
important physical features of the level instrument process. 
This model was simplified and used as a basis for a computer 
based steam generator level instrument. Although the model 
is based on fundamental properties and relationships, its 
accuracy and the limits of its application cannot be 
established fully until comparisons can be made with 


experimental data. 


1.1 Purpose of Simulation 


The optimum method for evaluation of the proposed steam 
generator level instrument would be to compare the outputs 
of a conventional level measurement device and the 
measurements of a microprocessor-based level instrument 
given identical inputs from an actual steam generator. 
However, since the concept is still being developed, the 
more prudent approach is to perform extensive lab tests 
prior to implementation of this smart instrument on an 


actual power plant. To evaluate the potential effectiveness 


22 





of the proposed smart level instrument, typical steam 
generator level pressure and power relations are established 
to provide simulated "measured" data for the physical input 
variables that the smart level instrument would require and 
for a simulated "known" level signal for the basis of 
comparison. While it is somewhat artificial to base 
measurement results upon simulated input data, this 
procedure does provide a qualitative and quantitative sense 
for the capabilities that a smart instrument would have and 


serves as a first step for future research. 


1.2 Smart Instrumentation 


In the context of this thesis, a smart instrument is one 
that includes a detailed model of the physical process being 
measured to enhance both the reliability and accuracy of the 
instrument. In addition, the model can give the instrument 
the capability to infer quantities that are not actually 
measured. For example, a temperature sensor could be used 
to control a heating coil for bringing a solution in a 
beaker to a desired temperature. A smart temperature sensor 
could detect the absence of the solution from the beaker and 
protect the heating coil by noting whether the rate of the 
temperature increase just after energizing the heater was 
consistent with that expected when the solution was present, 


and de-energize the heating coil appropriately. 


1.3 Background 


Accurate indication and measurement of plant variables are 
required for safe operation of nuclear reactors. 
Information recorded by instruments is used to advise plant 
supervisors about the operability and reliable operation of 


nearly all plant components, to provide inputs to protection 
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Systems, and to furnish input to automatic control systems. 
In the past, instrumentation was not sophisticated enough to 
enhance accuracy by using information about the process 
being modeled; often attempts were made to design 
instrumentation resistant to sources of signal degradation, 
but the design was usually based on worst case or very 
simplified assumptions about sources of measurement error. 
In addition, conservative estimates are generally made to 
account for the tendency of the transducer and associated 
circuitry to exhibit characteristics varying with time so 
that regular calibrations could be done prior to receiving 
erroneous readings. Recent attempts to integrate 
microcomputers into plant instrumentation and to employ 
signal validation techniques have illustrated the potential 
benefits of more sophisticated instruments (references S1 
and M2). Signal validation increases information 
reliability by combining information from redundant sensors 
and a detailed knowledge of the dynamics of plant systems. 
Detection of faulty sensor operation, a highly desirable 
function, 1S an important component of signal validation 
that can be accomplished by high speed comparison of 
redundant sensors in conjunction with knowledge of prior 


sensor data. 


The premise of this research is that microcomputers, when 
included as integral components of plant sensors, can 
perform the functions of signal measurement, validation, and 
fault detection to enhance instrument reliability and 
accuracy. Using microcomputers within each sensor enables 
one to program into the instrument an accurate model of the 
physically pertinent variables and their relationships which 
can enhance the accuracy of the instrument as well as 
provide a localized determination of the quality of data 
presented, also known as distributed processing. This 


information can be used to alert the operator to the onset 
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of conditions that could produce erroneous readings as well 
as extending the time between required instrument 
calibrations. In addition, redundancy and consistency 
checking of measured data can be done locally by a group of 


smart instruments. 


1.4 Motivation 


Accurate steam generator water level information is 
essential for safe operation of pressurized water reactors 
(PWR). The steam generators are the source of the steam 
used to drive the power turbines and are the heat sink for 
the primary coolant that flows through the reactor. 
Maintaining steam generator water level below the upper 
limit of the designated operating band is necessary to 
prevent carryover of liquid droplets to the turbines, which 
could cause severe turbine blade erosion and possibly 
tremendous damage to the turbine itself. Steam generator 
water level must not be allowed to drop below the lower 
limit of the designed band since this could uncover the 
generator U-tubes, resulting in inadequate heat removal from 
the reactor and significant thermal stresses to the U-tubes 
themselves, possibly leading to their premature failure. 
Clearly, it is desirable to enhance the reliability of the 
level information provided to plant operators and to the 
automatic control systems responsible for maintaining 
generator water level within normal bounds. Using a 
computer to record plant data is not a new concept, but the 
smart instrument would be able to interpret the measured 
data as well and by placing the computing resource as close 
as possible to the variables being measured, higher level 
computers utilized in the plant can be freed to perform 
integrating roles on the basis of input information of 


improved quality. 
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Chapter 2 


Steam Generator Description 


Each primary loop in a Pressurized Water Reactor contains a 
Steam generator, important and very large (typically 20 
meters high, 4 meters in diameter) components which are used 
to transfer heat from a primary coolant, flowing through the 
reactor core at high pressure, to a secondary coolant, which 
does not flow through the core, to produce the steam 


necessary to drive turbo generators. 


2.1 Steam Generator Internal Description 


There are three basic steam generator designs: vertical 
U-tube, once-through, and the horizontal U-tube. Since the 
vertical U-tube generator (Figure 1) is the most widely 
used, only its operation and level sensing configuration 
will be discussed here. The steam generator is a heat 
exchanger that enables the primary coolant to boil the 
secondary coolant without allowing the two fluids to mix. 
The hot, primary coolant just exited from the reactor vessel 
enters the steam generator at the inlet plenum, located at 
the lower end of the steam generator, passes upward then 
downward through the U-tubes, enters the exit plenum and 
returns to the cold leg reactor piping. Mixing of the inlet 
and outlet plenums is prevented by an isolation plate 
separating the two. Condensed steam, known as condensate or 
feedwater, enters the steam generator via an annular piping 
arrangement inside the generator called the feedwater inlet 


ring. Cold feedwater is prevented from impinging upon the 
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tubes by a tube bundle wrapper which forces the feed water 
through the downcomer region, so the feed water can be 
pre-heated by the hot water falling from the moisture 
separators. The feedwater flows up through the U-tube 
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Figure 1: Steam Generator Internals 


region where it is heated and transformed to steam because 


of the lower pressure of the secondary loop. The two phase 
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‘mixture passes upward from the U-tube region through the 
moisture separators, which allow the steam to pass but cause 
virtually all of the entrained liquid to fall back within 
the steam generator and mix with incoming feedwater, as 


mentioned above. 
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Figure 2: Steam Generator Level Sensing System 


2.2 Level Sensing 


Figure 2 shows a schematic representation of a typical steam 
generator level remote sensing system. It is analogous to 
the method used to measure core water level in a boiling 
water reactor (BWR) and the method used to measure the 


pressurizer level in a pressurized water reactor (PWR) 
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plant. The system depicted is sometimes called a cold 
reference leg system because the reference leg piping is not 
heated or insulated so it assumes the temperature of the 
ambient air in the reactor building. The downcomer level, 
the level in the area between the tube sheet wrapper and the 
generator vessel, is the variable actually measured and is 
considered the generator water level. That level is more 
easily measured and the two phase boiling process in the 
tube bundle region makes the level there ill-defined. As 
shown in Figure 2, an upper tap in the steam space is 
connected to a condensing pot and a reference leg. A lower 
tap in the liquid space is connected to a variable leg. The 
pressure difference, AP, between the reference and variable 
legs is measured by a differential pressure cell. The 
liquid level L may be calculated under the following 


assumptions: 


(1) the fluid densities in the reference and variable 

legs are known, 

(2) the liquid and vapor are separated, with liquid at 
the bottom and vapor at the top, 

(3) the reference leg is completely filled with liquid 
to a known height, 

(4) the density of the reference leg is constant, and 

(5) there are no velocity dependent contributions to 


the pressure differences. 


Figure 3 shows a schematic of the differential pressure cell 
arrangement, simplified to show just the details associated 
with the variable and reference legs. Based on figure 3, the 


differential pressure is approximately: 


5 








Figure 3: Differential Pressure Cell 


Ego -AD-upptcpeHg-n.Lg-p,g(H-L) (1) 


and so the level is: 


pa Pe Peg) -AP | IL (2) 
) 0 من ۵ 3 (۵< یی‎ ۵ 


where: De ۵ pressure in condensing pot; 
Dia â pressure in steam generator; 
p. à density of fluid in reference leg; 
Pac à density of fluid in the downcomer; 
p. 4 density of vapor in the steam generator; 
g 4 acceleration due to gravity; 
H 4 height of water in the reference leg; 
and 
L 4 downcomer level. 


Note that the level calculated by equation (2) is 
approximate since the height difference between the level of 


liquid in the condensing pot and the reference leg tap on 
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the steam generator is neglected. Otherwise, the numerator 
of the second term on the right hand side of equation (2) 


would be (p,H-p,H,.) where H, is the vertical distance 


between the reference and variable leg taps on the side of 
the steam generator. Since f, is so much larger than f,, 


this approximation introduces little error. 


2.3 Limitations of Conventional Level Sensing 


Generally, the assumptions used for level measurements and 
noted above are valid and equation (1) provides a good 
estimate of steam generator level during normal operations. 
However, severe transients involving rapid depressurizations 
such as a generator steam leak could lower the pressure 
above the reference leg enough to cause flashing in the 
reference leg, violating assumption 3. Variations in the 
temperature of the reactor building will affect the rate of 
condensation in the condensing pot and the average density 
of the liquid in the reference leg. Generally, reactor 
building temperature is moderately related to time at power 
the power level of the reactor, and the outdoor temperature. 
In addition, a casualty such as a steam leak from the 
generator or a primary coolant leak could alter reactor 
building temperature enough to make assumption 4 invalid. 
The diameter of the piping chosen for the condensing pot 
drain line is also an important consideration. If the 
piping is not large enough to accommodate the flow of 
condensing liquid from the condensing pot during worst case 
conditions, the pipe may become filled with liquid, 
resulting in a decrease in effective reference leg height as 
the water level in the condensing pot rises above the drain 
line piping connection to the condensing pot. This 
condition violates assumption 4. There have been 


well-studied examples of this phenomenon actually occurring. 
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۲:۲ ۳, the steam that flows from the steam generator to 
the condensing pot and condenses does so because the 
condensing pot is uninsulated, but it would not be accurate 
to describe the condensation as taking place at reactor 
building ambient temperature. Contact between the inner 
surface of the condensing pot material and the steam from 
the steam generator no doubt elevates the average 
temperature of the condensing pot surface above reactor 
building ambient temperature. Still, the portions of the 
reference leg piping in the vicinity of the differential 
pressure cell are most nearly at ambient for the reactor 
building. Clearly, assumption 4 is not really accurate and 
the density of the reference leg varies along the length of 


the reference leg piping. 
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Chapter 3 


Smart Instruments 


A smart instrument is an electronic device, possessing 
memory and controlled by a microcomputer, used to measure 
some plant variable. Its computational and decision making 
capabilities, the salient features distinguishing smart from 
conventional instruments, enable the designer to integrate 
within the instrument a model of the physical process 
affecting the parameter being measured. To implement a 
smart instrument, the designer must choose an appropriate 
model and consider the roles of measurement consistency, 
data redundancy, and filtering. Smart instruments have many 
potential advantages over conventional instrumentation, but 


there are some tradeoffs to be considered as well 
Translation 
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Figure 4: Smart Instrument Block Diagram 


3.1 Model 


Figure 4 shows a block diagram of a typical smart 
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Enstrument. It takes its inputs in the form of physical 
variables from some actual process and translates these 
natural phenomena into electrical variables like voltage or 
current. Conventional instruments process a transducer 
output signal in accordance with pre-determined assumptions 
selected by the designer. The processed transducer data is 
routed to a meter face and possible also connected to a data 
recording computer. While this has been a generally 
Satisfactory method for presenting plant variable 
information to operators, it has two drawbacks. First, the 
burden of evaluating the consistency and accuracy of the 
data presented is entirely upon the operator. This context 
evaluation burden is not onerous when the total number of 
variables that the operator has to monitor are small, but it 
becomes considerably more difficult in the control room of a 
commercial power reactor having thousands of instruments and 


potentially hundreds of alarms to monitor. 


Detecting instrument or transducer failure is especially 
complicated when the instrument continues to display a 
normal reading although its ability to respond at the 
extremes of its range might be severely degraded. The 
Second limitation of conventional instruments is that they 
are only accurate as long as the design assumptions remain 
valid. The microcomputer controlling the smart instrument 
enables the designer to program into the instrument an 
"understanding" of the numerical characteristics of the 
System being measured. The complexity and accuracy of the 
model chosen for inclusion into the smart instrument is an 
important consideration. The primary benefit of a good 
model is that it enables the instrument itself to monitor 
measurement results and alert the operator if the data 
received from the sensors is not consistent with either: 
redundant transducers measuring the same variables or with 


other sensor information that can be fed to the smart 
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instrument.  Judicious choice of the parameters to be 
monitored by the instrument and included into the 
computational model is essential. Selecting too few 
parameters and an oversimplified model risks realizing only 
a marginal improvement over existing instrumentation while 
an exceptionally complicated model that requires numerous 
interconnections with other sensors would require an 
enormously powerful computation engine that might dwarf the 
size of the original instrument, posing considerable 


debugging and robustness challenges. 


3.2 Data Collection 


It is important to establish some confidence that the data 
collected by the smart instrument is reliable. The 
important issues are: consistency of the detected signals 
with previous Signals from the same transducer, comparison 
of the detected signals with current redundant signals, and 
filtering the detected signal to reduce the influence of 


electronic or process noise. 


B 2.1 Consistency 


It is desirable that the smart instrument be capable of 
checking each measured data value for consistency with other 
variables. Consistency checking of instrument results is a 
key supervisory task performed by the plant operators and it 
is plausible that much data supervision could be performed 
by computer controlled instruments. The instrument should 
have memory of past measurement results so it can compare 
the current measurement with recent values of the same 
measurement or some average of the recent values. The smart 
instrument can check the current signal trend with the trend 


exhibited by other variables. The smart level instrument 
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could be programmed to regard a deviation of the current 
Steam generator level datum from the average of prior 
readings in excess of some threshold as normal if steam 
generator pressure is changing at some rate and abnormal if 
it is not. Incorporating a process model into the level 
instrument gives it the ability to determine if a particular 
measurement result is consistent with other plant variables. 
Since a smart instrument is envisioned as an aid to the 
operator, it should not require additional alarms or require 
Significant retraining to use (reference S1). The results 
of the consistency checks it performs could be presented to 
the operator by means of a small warning message on the 
instrument meter face or cathode ray tube display as well as 


being fed to a plant data logging computer. 


3.2.2 Redundancy 


Redundancy must be designed into any instrument that 
monitors an important reactor plant variable (reference M2). 
It would be a serious shortcoming for an instrument 
providing display of and control system inputs based on a 
plant variable measured by a single sensor since it could 
lead to undesirable operator or control system response 
either during normal operations (e.g. failure to act ona 
malfunction because it is not indicated by the instrument) 
Or in casualty situations (e.g. masking the casualty or 


giving inconsistent output). 


The smart level instrument should have at least 2 
independent level sensing signals as inputs. It is already 
a common procedure in power plant design to have 4 
independent electric level sensing channels per steam 
generator so this is not an overly restrictive requirement. 
The smart level instrument would record and compare the two 


level signals at each sampling interval. Level data points 
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could then be averaged over a suitable number (selected to 
establish a desired confidence interval) of sampling 
intervals. In addition, the smart instrument could be 
programmed with an allowable deviation from the average of 
prior signals, either from prior knowledge or using 
information recorded in real time, so that sensor data 
points that deviated in excess of the allowable amount from 
would not be included in the average. If a certain number 
of consecutive signal samples are rejected or some threshold 
of total signal samples over a certain time interval is 
exceeded, the operator could be alerted so the affected 


instrument channel could be checked. 


ES Filtering 


The instrument designer must anticipate that the sampled 
data presented to the smart instrument will be subject to 
noise degradation. The noise would result from interference 
from other electrical signals, process noise, and 


quantization noise. 


Radiation shielding considerations for the reactor and 
reactor building dictate that there be a minimum of 
penetrations in the secondary reactor shield to enhance 
shield structural integrity and reduce the possibility of 
radiation streaming. Minimizing the number of shield 
penetrations necessitates routing multiple electrical signal 
cables in bundles through a limited number of shield 
penetrations. If the smart instrument is located outside 
the secondary shield, it is a virtual certainty that cables 
connecting the smart instrument to its transducers will be 
placed in proximity to higher voltage and or frequency 
Signal cables for other systems so that there will be 
crosstalk from these other signals that will reach the smart 


instrument input. There will probably also be a significant 
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noise component at 60 Hertz due to power line interference. 


Process noise is composed of random signal variations 
introduced at the transducer input. The random component of 
the signal is produced by minor local variations in the 
measured values and the statistical nature of the physical 
interactions between the measured variables and the 


transducer. In addition, the smart instrument can be 





| programmed to evaluate sensor performance by interpreting 


the character of the process noise. The instrument could 
either be pre-programmed with maximum and minimum tolerable 
values for noise variations or could determine these limits 
through a weighted average of prior data. This expected 
noise deviation data could be used to alert the operator to 
an abnormal sensor noise signature and thus detect otherwise 


undetectable sensor failures. 


Quantization noise results from the inability of a digital 
computer to accept data in the continuous, analog form 
produced by the transducer. Before the smart instrument can 
process the sampled signal, it must first be converted to a 
digital representation. The resolution or accuracy of the 
analog to digital conversion is related to the number of 
binary digits, bits, used to represent the digital data. 
While quantization noise can be reduced by using an analog 
to digital (A/D) converter with more bits (at the cost of 
storing and transmitting more data), it is an important 
consideration for the designer, especially since the digital 
Signal values will be used in various computations of model 


parameters. 


To compensate for process noise and interference noise 
degradation of the sampled data, the smart instrument could 
be programmed to filter the incoming data digitally prior to 


processing in the model. It could also accept and record 
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the raw, unfiltered, digital data so it could be compared 
with the filtered data if necessary for conducting 
performance diagnostics. One could select a filter to 
remove a range of frequencies or reject a specific frequency 
where the majority of the noise power might be concentrated. 
A moving average filter would be useful for getting a good 


estimate of the mean signal value in the presence of noise. 


One other benefit of using a digital filter is that the 
filter characteristics could be easily changed without 
requiring a redesign of the instrument or the filter could 
be adaptive, possibly using a matched filter, to adjust the 
filter parameters in real time as the information is being 


received. 


3.3 Capabilities and Limitations 


A smart instrument’s capabilities include: 


‘Enhanced Reliability: a smart instrument can provide 
more reliable information than conventional sensors because 
it can include some extra variables in its model that enable 
it to provide reliable indication over a broader range of 
Situations, needs fewer assumptions, and can correct 
transducer drift if the designer has the ability to model 
Et. 


Calibration Reduction: smart instruments provide a 
real time means for objectively assessing the quality and 
accuracy of data presented, thus replacing strictly calendar 


based calibration requirements. 


-Variable Inference: in some situations, a smart 
instrument can infer plant variables that are not directly 


measurable, but are desired for increasing model accuracy, 
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EESessing instrument reliability, or for providing extra 


operational information. 


«Self Diagnosis: a smart instrument can be programmed 
with built in test (BIT) routines to perform diagnostics on 
its major components and alert the operator to specific 


failures, thus reducing repair time. 


«Sensor Communication: smart sensors can query other 
sensors for information on a common communication bus 
without requiring dedicated instrument to instrument 
connections. A control circuit could periodically poll 
sensors for data. Sensors could record the time when a 
certain reading or alarm occurs to be used for fault 
isolation. For example, pressurized water power plants 
usually have salinity cells at various locations in the 
condensate stream of the secondary system to detect the 
introduction of impurities into the condensate. 
Unfortunately, a large sea water leak upstream of several 
salinity cells will usually cause all the downstream cells 
to alarm nearly simultaneously so it is difficult to specify 
rapidly the source of the contamination, thus requiring more 
time to find and correct the source of the casualty. If the 
salinity sensors recorded the times when their alarms 
occurred, the operators might be able to identify more 
easily the source of the contamination based on the recorded 
alarm times. In addition, a communications channel between 
instruments could enable a control unit or protection 
circuit to query the instruments that provide it with input 
data, either to perform continuity checks or brief on-line 
diagnostics. Depending on the sampling interval and the 
hardware used, there could be sufficient time between data 


samples to do these checks. 


The limitations of smart instruments include: 
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-Model: the performance benefits discussed above are 
limited by the quality and accuracy of the model used by the 
smart instrument. The model must be chosen carefully so 
that it is robust and not subject to oscillatory or unstable 
E or. 

:Model Complexity: the designer does not have 
unlimited flexibility when selecting model complexity and 
variables to include. The more complicated the model, the 
faster the central processing unit needed and the more 
memory that will be required. In addition, robustness tends 


to be a more serious issue for more complicated models. 


‘Fault Tolerance: the smart instrument must be capable 
of withstanding power failures and a potentially harsh 
operating environment. Its data presentation algorithm and 
model implementation must be capable of recovering from 
unexpectedly large and small inputs without suffering 
failure. The designer might decide to use some kind of 
battery backup for the instrument, especially if it does 


| Some of its processing based on real time variable history. 


‘Interface: computers and computer controlled devices 
generally exhibit the tendency to produce and display data 
that has the appearance of truth even when receiving faulty 
inputs. It is important that the data be presented to the 
user in a form that enables him to make a determination of 


the quality of the data. 
Conclusions 


Because it is controlled by a machine that is capable of 
making complicated decisions based on input data and has the 


ability to adapt its performance to the measured data during 
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operation, the smart instrument is capable of considerably 
more sophisticated data analysis and manipulation than a 


conventional instrument. 
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Chapter 4 


Dynamics of a Steam Generator Level 


Sensor 


The model is developed in this research to relate the 
pressure in the condensing pot to that in the steam 
generator. The model permits calculation of the superficial 
vapor and liquid velocities and the liquid height in the 
drain pipe for the two phase flow from the condensing pot to 
the steam generator. None of these quantities is normally 
measured, but a level instrument that could accurately 
estimate their value dynamically could alert the operator to 
the onset of conditions that could lead to inaccurate level 


information. 


4.1 System Description 


The purpose of the condensing pot is to ensure the reference 
leg of the level instrument is always filled to a constant 
level. During normal operation, steam flows from the 
generator steam space to the condensing pot via a piping 
connection, called the drain line or tube in this paper. In 
the condensing pot, the steam condenses because the 
condensing pot is not insulated (as noted in earlier 
descriptions). Some of the condensed liquid keeps the 
reference leg full, but the remainder spills over from the 
condensing pot and into the drain pipe, returning to the 
Steam generator. The condensing pot is modeled as cylinder 
with its axis horizontal having the drain pipe connected to 


one of the circular bases. 
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A design aim for the drain pipe is that it be large enough 
to accommodate the counter-current vapor-liquid flow and 
remain smooth-stratified or wavy-stratified under most 
normal conditions. However, some combinations of drain pipe 
diameters and phase flow rates in the drain pipe can lead to 
a flow transition from wavy-stratified to either 
intermittent, in which waves at the surface are picked up to 
form a frothy slug which is propagated along the channel, or 
annular-dispersed, the formation of a vapor core with a 


liquid film around the periphery of the pipe. 


Rearranging equation (1) produces: 


TE Up ام‎ PAP ۵ م۱۳ ت91 لي‎ ۵ GE) 


As noted earlier, erroneous level indication could result 
from failure of the reference leg height to remain constant. 
As shown in figure 5, if the height of liquid in the 
reference leg rises above the design level (which would most 
likely a result of flooding in the drain pipe), the 
differential pressure transmitter would sense a decreased 
differential pressure since the effective H is actually 
lower. Since H is assumed by the instrument designer to 
remain constant under all circumstances, the instrument 
would interpret this change as a decrease in level. The 
effect of this condition on the instrument would be the 
presentation to the operator and to the level control system 


of an indicated level inconsistent with the actual level. 
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Figure 6: Condensing Pot Control Volume 


For the system under consideration, the steam generator, 
condensing pot, and the drain pipe, two control volumes are 


identified. The first consists of the condensing pot. The 





second is the drain pipe itself (between the condensing pot 
and up to the steam generator tap connection). They are 


depicted schematically in figures 6 and 7. 
The assumptions of the following calculations are: 


1. Each control volume is assumed to consist of saturated 
liquid and vapor in thermodynamic equilibrium under all 
operating conditions. 

2. The drain piping, usually consisting of various sections 
of horizontal and vertical piping runs, is simplified to be 
just one section of piping, inclined from the horizontal by 
an angle 8. This can vary from installation to installation. 
3. No significant concentration of non-condensible gases is 


present in the steam generator. 
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Figure 7: Drain Line Control Volume 


Thefmodel inputs are: p,, (pressure in the steam generator 
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at the tap adjacent to the drain pipe connection), p, 
(pressure in the condensing pot), a, (the average vapor 
fraction in the condensing pot), and h, (the liquid height 
in the drain pipe). The model requires the initial values 
fOr Peg, Per %, and h,. The inclination of the drain pipe 


from the horizontal, 0, is assumed known. 


The key difficulty when trying to incorporate features of 
the drain pipe flow into the overall model is how to treat 
and calculate the transient liquid height in the tube when 
steam generator pressure changes. For the purposes of this 
model, the connection between the two is horizontal or 
nearly so and the liquid surface in the condensing pot is 
approximately a free surface. It therefore seems reasonable 
to conclude that the characteristics of the flow at the 
transition between the condensing pot and the drain pipe are 
best described by open channel flow theory. The liquid 
flow input is varied during transients in the sense that the 
liquid height in the pipe at the junction with the 
condensing pot is the result of weir flow. The liquid mass 
balance in the condensing pot is a result of weir flow and 
the instantaneous rate of condensation in the condensing 
pot. The liquid height in the drain pipe closer to the 
Steam generator is a function of the liquid and vapor 


velocities and the pressure gradient along the pipe. 


4.2 Conservation Equations 


A first equation is obtained that defines conservation of 


mass in the condensing pot control volume. 





dM. ١ 
= A-m 4 
dt Dv Jv w ( ) 
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where M, ü mass of water in the condensing pot; 


p. à density of water vapor in the condensing 
pot; 
Jy â superficial vapor velocity of the vapor 


entering the condensing pot from the 
steam generator; 


A à the cross sectional area of the drain 
pipe; and 
um á the mass flow rate of liquid leaving the 


condensing pot in the weir. 


Note that the superficial liquid velocity in the weir is not 
necessarily the same as j,, the superficial liquid velocity 
in the drain pipe close to the steam generator, at least 
during transient conditions when the pressure of the steam 


generator and the condensing pot are changing with time. 


Conservation of energy, in terms of enthalpy, in the 


condensing pot control volume gives: 


dH d ۱ 
DE -vE p, j AH h, HÀ (5) 
where H: ۵ the total mass-enthalpy product in the 
condensing pot; 
M. å the total volume in the condensing pot; 
A à the total cross sectional area of the 
drain pipe; 
ES à the enthalpy of the vapor entering the 
condensing pot from the steam generator; 
H, â the enthalpy of the liquid leaving the 
condensing pot; and 
Ò â the heat loss from the condensing pot to 
ambient. 


The total mass and enthalpy in the condensing pot are: 


M.-[p,(1-«,) * p,a,] V, (6) 
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7 OTO O HCL (7) 


where Pe ۵ density of the liquid leaving the 
condensing pot. 


In the drain line control volume, conservation of mass 





gives: 
aM. 
=n, + Pp ÍA (8) 
de wT Pete 
where M, à total mass of water in the drain line or 


tube and the influence of the vapor 
phase in the pipe has on the change in 
mass in the pipe has been neglected. 


As noted earlier, it is assumed that the transient flow 
condition at the entrance to the drain line, at the 
interface to between the condensing pot and the drain line, 
can best be modeled as a weir. From open channel flow 
theory, reference Cl and Ll, for steep channels, the flow is 
critical at the channel inlet, thus the flow energy is 


minimized there, yielding the critical flow condition: 


wie 


8 
aa.) مایت کت‎ 
a, 


=1 (9) 


Bo that: 
; gD,cosO 
10, 70,7,,A70,(1-4,) A u u (10) 
t 
where Blais à superficial liquid velocity at the weir, 


0 â vapor fraction at the boundary between 
the condensing pot and the top of the 
drain pipe, 


0 à angle of elevation of drain pipe from 
Horizontal, 
g â acceleration of gravity, 


39 





D. - diameter of the drain pipe, and 
aq 3 energy or coriolis coefficient. 


a, only applies in the vicinity of the weir and is not the 
same as the void fraction in the pipe, Q,, closer to the 
steam generator. The energy coefficient, a,, varies with 
the liquid depth and is usually determined empirically, is 
set equal to 2 for simplicity since, according to reference 
Cl, experimental values of a, are in the range of 1.03 to 
1.36 for fairly straight, prismatic channels. Since the 
precise degree to which the connection between the drain 
pipe and the condensing pot approximates an open channel is 
not known, the value for a, is taken to be 2 as a worst case 


estimate. 

Differentiating (3) and (4) 

dM, 1 m ape ( = da. 
Ze 


E lI grav de 


c 

















AM. _ Op, 06 apa au. 
E | ID. (1 & e) "ép. * Vee * (Py Pe) Ve (11) 
and 
ar OH dp. , OH.| d«, 
dt | 9. dt |9%,) dt 
Eee pH) (1-a eh) a lV Pe + (p H,- pH) V Me (12) 
at OD. c Op. c| € dt vv t ce de 
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۴۳۳ ۰ ۳۳1۳9 (11) and (12) into (4) and (5) 





Op, dp, dp. da, 1 
O A Ve gE PAi (13) 





da ۱ ; 
۲ (p, H, - p, Hy) Ve-dp C8vJVyAH, - fu H,- Q (14) 


The total mass in the drain line, M, is: 
M.-A.L.p, (15) 


where L. ۵ the total length of the drain pipe; and 
A,=(1-a@,) A. 


Relate M,, the total mass of water in the drain line, to h,, 
the liquid height in the drain line, through the following 


geometric relationships (see figure 8): 


4 ومع‎ | 1-2 2 (16) 
t 


2 
E ($,-cosó,sin$,) (17) 


Exbststuting (16) into (15): 


2 
«AA ($ - 054 ,6 11 ( (18) 
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Figure 8: Relation Between h, and > in the Drain Pipe 





differentiating (14) 














d 
+ sin*$, E + asin, St - 2 52124 م‎ ۳ (19) 
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recall 











d 2 p 
— (cos x)-- 
dx /1-x? 
Using the above formula to differentiate (16) produces: 
d dh 
X : dc 1 
fy h, 
D4|4—|1-— 
pe D 
EN 1 ۱۱ (21) into (20) 
sin cos"(1-2 2% 
eM, De, Dc (22) 
dt 2 ۳ Eu 
a a 
ee 
D. De 


Eubstrtuting (22) into (8) 


h 
sin?^cos !|1-2— 
D, 


Eombining (13), (14), and (23) into a matrix equation: 


—L.D, 


dh, , 
E سس‎ "mhp + Pj (23) 


ele 


Oli 
1 
| b> 


ct 


e 
—_x=f 24 
At = 
where: 
A. A5 0 
A=|42, 422 0 (25) 
0 0 A, 
and 
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Po p AU 
x= “| (26)  f-p,j,AH,-m,H,-Q )27( 
h, Ih,* 9,7,À 


E Eusrres of A are: 
| 








a, x 








A227 (pH, -P A) Ve 


1 -1 A, 
D SIID^cos us 

= t 
DES Lp, - 


a 
t D, 


1 


4.3 The Functional Dependence of f 


EE 1 jg. and 7, 


While equation (27) identifies the components of f used in 


equation (24), it is useful for computation to clearly 


Au^ 
ar a V 
A= I (oH) (1-« o E 0,۷ 
1 OD. $t € Op, ۷ € Cc 


establish the functional dependence of the individual 


components of the f. The method for finding the 
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superficial velocities, j, and jı, is adapted from the 
approach used by Taitel and Dukler, reference T1. Figure 9 
shows a smooth equilibrium stratified flow on which the 


following equations and geometric relationships are based. 


A momentum balance on the liquid and vapor phases yields: 


To Condensing Pot 


+x direction 





Figure 9: Two Phase Flow Geometry In Drain Pipe 


Ay 1, 8,71, S,-A p, gsin-o (28) 
A E <q کی‎ + 5-2 in8-0 (29) 
az 4 9 OATS 
where A, à cross sectional area of the vapor phase; 
A, ۵ cross sectional area of the liquid 
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phase; 

wall to vapor phase shear stress; 

wall to liquid phase shear stress; 
interfacial shear stress; 

wetted perimeter of vapor phase; 

wetted perimeter of liquid phase; 
interfacial wetted perimeter; and 

angle of pipe elevation from horizontal 
(see figure). 


b p لا‎ pC CED pp 


using the nomenclature from Tl and C2. 


From reference Tl, the geometric quantities are: 


A 
S¡=4D, where pecos 1 =2 3] 
t 

S,-D,sin$ 


S,=D.(n-6) 


A,=A-A, 


The hydraulic diameters are: 








4A 4A 
دد‎ and D,—3 
DE Sy 
the velocities are: 
A 
U - 22 and go I 
i v 


The liquid and vapor friction factors, using a Darcy 


Weisbach representation, are: 


where "n ae dynamic viscosity of vapor and liquid, 
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۳۵ 
۱ E .184; and 


n = 2 


v,é 


for smooth stratified flow, 


f,=f 


| مک‎ 
The shear stresses are evaluated from: 


| 
عبر‎ eee |g = Fv Petey 


wi 8 8 


۱ هد‎ £p ,|u,- ui (u,- uj) 
j ti 
8 
Equations (28) and (29) are 2 non-linear equations in 4 


unknowns:  j,, ji, hp, and P’ (where P =£). By defining J, 


and ją as functions of h, and P', that is: 


pc Eng pr m) 


equations (28) and (29) can be iteratively solved for j, and 
J, in terms of h, and P'. Note that P’ is the difference 
between condensing pot and steam generator pressures divided 
by the length of the drain pipe. Equations (28) and (29) 
are iteratively solved for j, and j; using known values of h 
and P’. Note, the equations used to solve for the 
superficial velocities are in error because the shear 
stresses above were calculated using a factor of 2 in the 
denominator, not 8. This error was discovered too late to 
re-run the simulations. The resulting sets of data were 
tabulated and used to correlate the superficial phase 
velocities given known values of liquid height in the drain 
Pipe and pressure gradient along the pipe. The superificial 


velocity data is shown plotted in figures 10 and 11. 
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Superficial Velocities 
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Differential Pressure (Pa/m) 
Figure 9: Superficial Vapor Velocity Correlation 
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Figure 10: Superficial Vapor Velocity 
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Superficial Velocities 


w 
- 
~ 


Superficial Liquid Velocity (jl [m/s]) 





-600 -400 -200 0 200 400 600 
Differential Pressure (Pa/m) 
Figure 10: Superficial Liquid Velocity Correlation 


Figure 11:  Superficial Liquid Velocity 
4.3.2 The Weir Mass Flow Rate m 


m 


wW 


y 18 known from equation (10) using the critical flow 


condition, (9), but the relationship between œ, and QA, the 
vapor fractions at the condensing pot-drain line interface 
and in the condensing pot, respectively, must be more 
clearly identified. 

Note: (04 is the vapor fraction at the boundary between the 
condensing pot and the top of the drain pipe, it is not the 


same as the vapor fraction throughout the drain pipe. 


The weir mass flow rate is: 


1,7 j 1, A0, (10) 
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From the definition of Q and figure 8, the following 
equations relate h, to 0: 


A 
(مع-1)‎ (31) ۱ 
2 Note that these equations 
apply in the drain line and in 
and the condensing pot 


2-2) (6-cos@ sing) (17) 


Rearranging (16) 


Le =1-cosd, (32) 


c 


and solving for h,, (water height in the condensing pot) 


h,,7 5 (1-cos$,) (33) 


ulm 


substituting (17) into (31), and letting D-D,, 06=0,, and aza, 
yields 


(1-0) = (2) [decsimbecoóho) (34 


so that: 





$ —Ssing$,cosp,= (1-4.) Co =(1-a,) (35) 


c 


2 


If A. is known, (35) is used to find ọọ., then (33) is used 
to find hie. Knowing hier h, (liquid height at the weir, see 
figure 12) is found by: 
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If >یرط‎ then h,=0 (36) 





If ho — = then  h,-D, (37) 











Ww 


then h,=h,- | um (38) 





Figure 12: Relation Between h,. and h, 


Where the small angular dependence of h, is neglected. 
Clearly, only a certain range of h,, values will produce 
values of h, in the range from O0 to D,. If hh, is above the 
top of the drain pipe, the pipe entrance will always be 
filled. If h. is below the bottom of the drain pipe 


connection to the condensing pot, the pipe will be empty of 
liquid. 


Once h, is found, %, the vapor volume fraction at the weir, 
is computed from: 


Now that œ is known, Myr the liquid mass flow rate at the 
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bancos {1 Ze (39) 


t 


AN $,-sin$,cosó, 


x = (40) 


1-a,= 


weir present as the condensing pot-drain pipe interface, is 


found using equation (10): 


S (10) 


a 


1,7A.p, (1-0) 4 


4.3.3 The Heat Transfer Rate Q 


The heat lost from the condensing pot to the reactor 


building can be expressed as: 


Q=UA (Toen Tair) (e 


UA is calculable from the geometry of the condensing pot and 
some empirical relations for the film heat transfer 
coefficients of the liquid condensing on the inner surface 


of the condensing pot and of the air on the outer surface. 


For a cylindrically shaped condensing pot, 


A A ا‎ ASA (42) 
1 Tout 
1 en y 1 


1 سے‎ oM 
LN IL Kk. ARNT ERB 


SOUL er alr 


where r. â inner radius of the condensing pot; 
oe. ae outer radius of condensing pot; 
he ra film heat transfer coefficient for the 


water vapor condensing on the inside 
surface of the condensing pot; 
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hm A film heat transfer coefficient for air 
in the reactor building; 


re ۵ thermal conductivity of condensing pot; 
and 
B. à inner length of condensing pot. 


liquid air 
film condensing film 


Ir outside 
in condensing 
inside 
condensing 
pot 





Figure 13: Condensing Pot Wall Temperature Profile 


Figure 13 illustrates the typical temperature profile from 
the condensing pot to ambient. An empirical expression for 
Hoan, from reference M3 based on the Nusselt method for 
predicting the rate of heat transfer from a vapor condensing 
in a laminar film and running down a vertical surface, is 
used to get an order of magnitude estimate for the size of 


the film heat transfer coefficient of the condensing water 


vapor: 
4k Heo GL, _\ 4 
h ===VNu, where Nu, سس‎ ۳ | 4 43) 
stm 37, f t E i 
where: hs à evaporation increment; 
thermal conductivity of water; 
u ۵ dynamic viscosity of water; 
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[e à wetted length of condensing pot above 
normal water level; 


g à See ۵ رت‎ OF GEAVIEY; 

p à density of liquid water at condensing 
pot pressure; 

u G Temperature corresponding to saturation 
pressure in condensing pot; and 

ji à Temperature of inner surface of 


condensing pot; 


assuming: 


1. Vapor which flows parallel to the surface flows so slowly 
that its drag on the condensed film can be neglected. 

2. Motion of any liquid element in the condensed film is 
governed solely by a balance of the gravity and viscous 
forces acting on the element, hydrostatic pressure forces 


and inertia forces can be neglected. 


Despite the fact that the Nusselt correlation applies 
Specifically to vertical plates, it is useful for getting an 


estimate of the magnitude of ham- 


To estimate h,,,, a simplified expression for the heat 
transfer coefficient is used, obtained from reference M1, of 


the form: 


b 
neal 27) (44) 


S 


where A and b are constants, depending on geometry and flow 
conditions, and L, is the significant length, also depending 
on type of geometry and flow. Using the values suggested by 
reference M1 for an horizontal cylinder and laminar flow, 
A=.27, b=1/3, L, is the diameter of the condensing pot, and 
AT is T,, the outer surface temperature of the condensing 
pot, minus T,,, the temperature of the ambient air in the 
reactor building. Heat transfer by radiation from the 


condensing pot to the reactor building, a factor that is 
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less than 20$ of the overall heat transfer coefficient of 
Bar, ıs ignored. Including the contribution due to 
radiative heat transfer would slightly increase the value of 
Eheccalculated h,,. 


The difficulty in computing both the film heat transfer 
coefficients is that the value of each coefficient depends 
on the magnitude of the temperature difference across the 
respective films, which one cannot know until the heat 
transfer coefficient is known. The heat transfer rate from 
the condensing pot to the air of the reactor building, 


assuming that an estimate of j, is available to compute ih, 


can be found from: 


Q-i,hg, (45) 


where: Ih, à mass flow rate of liquid condensing in 
the condensate pot. 


If the thermal resistance of the condensing steam film is 
known, the temperature difference across the film can be 


expressed by: 





ando uu 1 


T,-T,)= AAA Bl 
(Term Tin) sa NTR 


Q 
Re 


stm 


where R: „ô thermal resistance of the condensing steam 


film, which is a function of the temperature 
difference across the film. 


Equation (46) can be solved iteratively to find the thermal 
resistance of the condensing steam film and thus h,,,. The 


Same procedure is used to find h,,. 


Figure 14 shows a typical condensing pot and drain pipe 


arrangement. Using a pressure of 7.34 MPa to compute hg, 
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and an estimate of .02 for j,, m, and Q can be computed 


using equations (10) and (45), respectively. The wetted 
length, L,, in equation (43) is taken as one half of 3.5 
znches, the length of the condensing pot shown in figure 14, 
since the condensing pot is nominally half filled with 


liquid at all times. 


3.5 inches 


2 inch dia 
Sched 80 condensing 
pot 


] inch 


reference leg 





Figure 14: Geometry Of Condensing Pot 


Taking k = .626 W/(m-K), g = 9.807 m/sec’, L, = 44.5 mm, r,, 
= 24.6 mm, ü = 90 uPa-s, p, = 729 kg/m?, an iterative 
solution of equation (46) and using equation (43) for hatm 
yields a value of h,, = 5400 W/(m’-K). An iterative 
solution of equation (46) using the expression for h,i: of 
equation (44) and assuming that h,,, is so large in 

that its effect of the total 


differential temperature is negligible, gives h,, - 12 


comparison with h,,, 
W/(m^-K). Returning to equation (42), using the above 


values and r,,-27.4 mm, k, = 50 W/(m-K), one finds: 


out 
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———————— 20.000663 
ARI: neum 
ii = 
out 
=0.008 
ZU, de : 
1 
سب‎ 21 00 9 
ZI E ir 


These results clearly show that the air film heat transfer 
coefficient dominates in the denominator of equation (42). 
For the purposes of simulation, a constant value of h,, is 
used because the actual value during transients at normal 
operating pressure, a function of the pressure in the 
condensing pot, will not vary too far from the nominal value 


calculated above. 


The resulting functional dependence of Q is 


OA BIT) (47) 


Having specified the functional forms of each component of 
f from equation (24), it is easier to see how to compute 


the necessary partial derivatives that will be needed to 
solve (24). 


4.4 Calculations 


The equations to be solved are of the form: 
Gina 
—X=f 
{me 
and 


3 





3 
ses, 


where: 


API A ASMA) 
E )را لدم از لصا‎ JA (a )H(p.) -OQ(pp) 


£,=M, (0 ¿) ^P; (D) IML, h) A 
and 


_ (Dc- Pag) 


p! 
Le 


Converting to a one-step, implicit difference of the 


linearized equations: 


A Es (21-22) 


At Ox 
and we seek xX"*1: 


-1 
az (az) f?^«X (48) 
Ox 





214.1 The Jacobian of the f vector 


The Jacobian of f is computed using numerical 


differentiation, where the Jacobian component in row i 


column j is: 
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0261 در‎ AX; 


Of course, it is important to choose the A value small 
enough to avoid errors due to nonlinearity of the subject 
functions and large enough to avoid roundoff in getting a 
meaningful average value for the derivative. Several plots 
of the functions f,, مرگ‎ and f, were computed to estimate 


adequate A values. The values used are: 


AD,730 ۵0,2۶ ۰0005 ۵2, 2 5 


The complete Jacobian, 2 of f is: 




















af, af, af, 
dp. da. An, 
af = Org REO. 
(32) | dp, ða, Oh, 
af, af, of, 
[7 da. Oh, 


For details on the software used to solve equation (41) for 
various initial conditions of X and pressure transients in 


the steam generator, please refer to the Appendix. 


4.4.2 Geometry and Initial Conditions 


For the purposes of these calculations, the following 


average parameters were used: 


Length of drain piping (L,) =2.5 m; 


Reactor building ambient air film heat transfer coefficient 
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(h.,.)=11 W/ (m’-K); 


Inner radius of condensing pot (r,,)=24.6 mm; 


air 


Outer radius of condensing pot (r,,) - 27.4 mm; 


Ambient temperature of Reactor Building (T,,)-50?C 


To calculate the initial or steady state values of X, a 
initial value for p, is assumed. Knowing this variable and 


the ambient reactor building temperature, equation (46) is 
used to find Q, the heat rejected by the condensing pot to 


ambient. With Q known, the liquid mass flow rate from the 
condensing pot the steam generator is found from equation 
(45). 


Equation (10) is used to compute the superficial phase 
velocities, presuming that the liquid and vapor mass flow 
rates are equal at steady state. Once j, and j, are known, 
equations (28) and (29) are solved iteratively until values 
for h and P’ are found. With a calculated P’ and 2 
EEUUSSOÍL Ob. p. can be computed and the process repeated 
until convergence is obtained for a value of p.. This is 
how the steady state values of X, composed of p,, œ. and بط‎ 
were computed. This process is diagrammed schematically 


below: 


P =0Q=m,=h _=0. 


la 


The above procedure for computing the initial conditions 
based on a combination of steam generator pressure, p,,, an 
assumed pressure gradient between the condensing pot and 


Steam generator, P', and the geometry of the condensing pot. 
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It is, however, possible to calculate a liquid mass flow 
rate from the condensing pot to the steam generator, My, 
that makes it impossible to solve the momentum balance 
equations for the liquid and vapor phases in the drain pipe, 
equations (28) and (29). Recall that at steady state, 
liquid superficial velocity in the drain pipe, j,, is 
related to m, strictly by the liquid density and the pipe 


cross sectional area (equation (10)). Reference W1 notes 
BEES plot of j, and j, for constant liquid height, Br, 
figure 23 for the tabular data calculated, enables one to 
identify graphically the regions of permissible operation in 
the drain pipe. Positive j, and positive j, represents 
cocurrent flow up from the steam generator to the condensing 
pot using the positive sense shown in figure 7. Negative j, 
and negative j, represents cocurrent down flow from the 
condensing pot to the steam generator. The region of 
positive vapor flow and negative liquid flow between the 
axes and the locus of points joining tangent lines to the 
lines of constant vapor liquid height, called the flooding 
line, represents the limited region of counter current flow, 
where solutions are possible (see figure 23). Operation 
above and to the left of the flooding line is not possible, 
because the drain pipe is already flooded at that point. 

If, while calculating the steady state solution, one picks 
assumed conditions of pressure gradient, P', that results in 
superficial velocities that are beyond the flooding line, 
one must continue to select a lower pressure gradient and 
recalculate the superficial velocities until physically 
realizable values of the superficial velocities are 


obtained. 


4.4.3 Model Limitations 


61 





The tabular data, values of j, and j, for known values of h, 
and P”, must be calculated before implementing the other 
model equations. The limits of the allowed liquid heights 
in the drain pipe are determined by the drain pipe diameter, 
but the range of potential pressure gradients cannot be 
known until after the model is implemented and some 
solutions calculated. The possibility exists for steam 
generator pressure changes to occur rapidly enough for the 
pressure gradient between condensing pot and steam generator 
encountered during transients to exceed the maximum or 
minimum values of P' used to compute the tables of j, and j, 
data. To account for this possibility, there are two 
options. First, perform extensive tests of the model using 
the worst case positive and negative steam generator 
pressure transients and expand the range of drain pipe 
pressure gradients used to generate j, and j, by trial and 
error. This is fairly easily, though tediously, 
accomplished and only increases the memory that the smart 
instrument must have. Second, if the smart instrument is 
capable of performing the calculations rapidly enough, the 
designer could elect to solve equations (28) and (29) at 
each time step rather than relying on tabular data computed 


in advance. 
2 4.4 Solutions 


The independent variable in the implicit difference equation 
(45) is p,,, steam generator pressure. Equation (48) is 
solved for the state variables condensing pot pressure, p,, 
condensing pot vapor volume fraction, @, and drain pipe 
liquid level, h,, for both constant steam generator pressure 


and up and down steam generator pressure ramps. 


4.4.4.1 Constant Pressure 
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Table 1 show the results of the solutions computed at a 


steam generator pressure of 7.34 MPa and various condensing 


pot and drain pipe diameters. 
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Steam Generator 


Pressure = 7.34 


Condensing Pot 


Diameter (mm) 


Drain Pipe 25.4 ZARA T9 27.4 9 
Diameter (mm) 


. 741 4.0 . 724 Z9 


Condensing Pot NU 3 
Vapor Volume 


Praetion e, 


| Drain Pipe Liquid 


Level h, (mm) 





Table 1: Steady State Solutions for Constant Steam 


Generator Pressure 


A nominal steam generator pressure of 7.34 MPa was chosen 
since many commercial steam generators operate between 7.6 
MPa and 6.9 MPa when steaming at normal operating pressure. 
The piping dimensions are varied to see the influence of 
dimensional tolerances on the stability of equation (48). 

As the data in table 1 shows, P' is the most sensitive state 
variable with respect to the diameter of the drain pipe. A 
7.9 percent increase in drain pipe diameter reduced the 
pressure drop by 21.9 percent, the vapor volume fraction 
rose by only 4 percent and the liquid level decreased by 
only 11.5 percent.  Physically, it makes sense that 
enlarging the drain pipe diameter would reduce the wall 
friction forces, thus reducing the pressure drop. The vapor 
volume fraction would also decrease since expanding the pipe 
diameter reduces the liquid level in the pipe. The last 


column of table 1 indicates that condensing pot vapor volume 
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fraction is the only state variable sensitive to condensing 
pot diameter and only weakly at that. One should also note 
that the pressure drop between the condensing pot and the 
steam generator is an extremely small fraction of steam 


Benerator pressure, about .0025 percent. 


4.4.4.2. Pressure Ramps 


Table 2 shows the steady state solutions reached for equal 

and opposite steam generator pressure ramps of 300 kPa, both 
Starting from 7.34 MPa. Pressure ramps were used as a first 
approximation for the generator pressure changes that result 


from changes in turbine loading. 


Down Ramp | Up Ramp 


= 300 4300 


(Pe Pag) ۱ -1 8 4 . 6 -179.42 


Condensing . .7144 RZ 


ESE Vapor 
| Volume 


۱۳۳ ۱۵1 O% 


Drain Pipe 
Liquid Level 
hy (mm) 





Table 2: Steady State Solutions for Pressure Ramps 


The trend of the steady state data corresponds well to what 
is predicted by the physical situation. For the negative 
pressure ramp, lower pressure in the steam generator results 
in less condensation in the condensing pot and reduced heat 


transfer from the condensing pot to ambient since the 
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condensing pot is assumed to be saturated and thus the 
temperature of the water vapor in the condensing pot is 
lower for the lower pressure. Reduced condensation in the 
condensing pot results in a lower liquid level there which 
causes the liquid level in the drain pipe to be lower (from 
equations (36) to (38)). The opposite trend is observed for 


the positive pressure ramp. 


As mentioned in section 4.5.3, since the values of P' used 
to solve equations (28) and (29) for j, and j, were 
specified in advance of computing the solution to equation 
(48), the possibility exits that transient values of P' 
outside the +600 Pa/m range those used might result. 
Preventing this form occurring constrains the maximum 
allowable rates of steam generator pressure change that my 
model could accept. As shown in table 3, the maximum 
allowed steam generator pressure gradient with respect to 


time is also different for up and down pressure transients. 


— = -— ——À— — ہے‎ pouxuu——————MÓÁÁ At m t — 


Generator Pressure Ramp Maximum Rate of Change 


Up .42 kPa/sec 





Table 3: Maximum Generator Pressure Gradients 


For the negative pressure gradient, figures 15 through 18 
show the transient steam generator pressure, the 
differential pressure between the condensing pot and steam 
generator, the vapor volume fraction, and liquid height in 
the drain pipe. Note that both the differential pressure 
and the drain pipe liquid height asymptotically approach 
constant values during the transient, but rapidly converge 


toward their steady state values after the transient. It is 
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particularly interesting to note that the these solutions to 
equation (41) show that the steam generator level model is 
capable of directly inferring the proximity to a flooding 
condition because an instantaneous value of liquid height is 
calculated during the transient. Other methods of 
predicting flooding rely on correlations between the 
Superficial liquid and vapor velocities in the pipe since 


liquid height in the pipe is not usually calculated. 


x106 -300 kPa Pressure Change In 3.75 Minutes 
T35 
ES 
725 
= 72 
e 
> 
2» 
E 7.15 
T 
7.05 
> 
-50 0 50 100 150 200 250 300 
Time (sec) 
Figure 15: à -300 kPa Pressure Ramp, p,, 
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Differential Pressure for -300 kPa Pressure Change In 3.75 Minutes 
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Figure 16: -300 kPa Pressure Ramp, p.-p,, 
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Figure 17: -300 kPa Pressure Ramp, œ 
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x10 Drain Line Liquid Level for -300 kPa Pssr Chg In 3.75 Min 
11 ۱ 


10 
9 
E s 
a 
7 
6 
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-50 0 50 100 150 200 250 300 
Time (sec) 
Figure 18: -300 kPa Pressure Ramp, h, 


Figures 19 through 22 show the transient steam generator 
pressure and state variables for the positive steam 
generator pressure ramp. As discussed, the positive 
pressure ramp had to occur over a time interval 3 times 
longer than that for the negative ramp in order to reach a 
solution. It is interesting to note that for the positive 
pressure transient, the transient liquid height in the tube 
increases just as it did for the negative pressure 


transient. 
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x106 +300 kPa Pressure Change In 13.33 Minutes 
T65 


T6 


7.55 


7.5 


Psg (Pa) 


7.45 


7.4 


7.35 


-100 0 100 200 300 400 500 600 700  Á 800 900 


Time (sec) 


Figure 19: +300 kPa Pressure Ramp, Pp,, 
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Differential Pressure for +300 kPa Pressure Change In 13.33 Minutes 
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Figure 20: +300 kPa Pressure Ramp, p.-p,, 
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Figure 21: +300 kPa Pressure Ramp, Q 
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۲10-3 Drain Line Liquid Level for +300 kPa Pssr Chg In 13.33 Min 
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Figure 22: +300 kPa Pressure Ramp, h, 


For the negative pressure ramp, the rising trend in drain 


pipe liquid height can be explained as follows: 


INPP.) changes from negative to positive as 

soon as steam generator pressure begins to decrease 
because condensing pot pressure cannot change as 
rapidly as steam generator pressure. 

2. To balance the momentum equations, (28) and (29), 
j, must become less than zero (see figure 10, j, must 
be less than zero when p,-p,, becomes greater than 
zero). 

3. Since the initial change in p, is small, QA, remains 
constant or changes slowly at the start of the 
transient, suggesting that h, is initially constant. 

4. Looking at figure 23, if j, becomes negative and h, 


remains constant, j, must become more negative as well, 


"m 





possibly by a factor of 2 to 3. 
5. Figure 11 shows that as j, becomes more negative 


for a fixed value of P', h, rises. 


For the positive pressure ramp, the rising trend in drain 


Pipe liquid height can be explained as follows: 


iP. Ps, becomes even more negative as soon as steam 
generator pressure begins to increase. 

2. To balance the momentum equations, (28) and (29), 
j, must become more positive (see figure 10, j, must 
increase when p.-p,, becomes more negative). 

3. Since the initial change in p, is small, «a, remains 
constant or changes slowly at the start of the 
transient, suggesting that h, is initially constant. 
4. Looking at figure 23, if j, becomes more positive 
and h, remains constant, j, must become less negative. 
5. Figure 11 shows that as j, becomes less negative 
for a fixed value of P', h, rises. h, would actually 


decrease if the P' only became slightly more negative. 
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One Dimensional Horizontal Flow-const h 
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Figure 23: Flow Map for Drain Pipe 


1.5 Flow Transitions 


A major goal of the steam generator water level model is to 
create an instrument sophisticated enough to alert the 
operator of the approach to flooding conditions in the drain 
pipe. If flooding does not occur, the pressure difference 
across the drain pipe remains very small so no correction is 
required. Recall that the drain pipe returns the overflow 
of condensed steam from the condensing pot to the steam 
generator and should be sized large enough to transport the 
maximum expected condensation flow without undergoing a flow 


transition from stratified flow to annular or intermittent 
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flow. Such a flow transition changes the height of the 


reference leg and produces erroneous level readings. 


4.5.1 Taitel and Dukler 





Taitel and Dukler (reference T2) developed a theoretical 
model for predicting the onset of flow regime transitions in 
near-horizontal piping that agrees well with experimentally 


observed results. 


The procedure was first to visualize a stratified liquid and 
then the mechanism describing the change from stratified 
flow. Beginning with stratified flow, as the liquid 
velocity increases, the liquid level rises and a wave is 
formed which grows rapidly, tending to block the flow. At 
lower vapor velocities, the wave crest forms a complete 
bridge and slug or plug flow results. At higher vapor 
velocities, there is insufficient liquid flowing to maintain 
or even form the bridge and the liquid in the wave is swept 
up around the pipe to form an annulus with some entrainment 


if the vapor velocity is high enough. 


Taitel and Dukler found that the criterion for the 
transition from stratified to intermittent and to annular 


dispersed flow is: 


SS 


!|»1 (48) 
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where F is the Froude number modified by the density ratio: 


F= Py E. m (49) 
| ۵ م۵‎ 8 


C2 is estimated as 1-h,/D, and the dimensionless values in 
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equation (48), signified by a -, are found as follows. 


À,- [cos (28) - (28-1) A - (28,-)?] 


A= [t cost (28) - (25-1) /1- (2&7 


2-22 
p.24 
^ fH 


Clearly, if one knows the liquid height in the drain pipe 
and the superficial vapor velocity, j,, the proximity to the 
Froude number characteristic of the flow transition can be 
calculated from equations (48) and (49). The chief obstacle 
to directly implementing a scheme to evaluate the Froude 
number is that the superficial velocities and liquid height 
in the drain pipe between the condensing pot and the steam 
generator are not measured. Despite the numerical error 
that was made when computing the shear stresses in equations 
(28) and (29), the results of section 4.4 suggest that it is 
theoretically possible to infer these variables from 
knowledge of the steam generator pressure and piping 


Geometry. 
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Ss Drift Flux Model 


As discussed, the flow in the drain pipe is normally two 
phase, counter current flow, with steam flowing from the 
steam space of the steam generator to the condensing pot and 
condensate returning from the condensing pot to the steam 
generator, illustrated in figure 9. The drift flux model, 
presented in reference Wl, a separated flow model focusing 
on the relative motions of the individual phases in the 
pipe, is well suited to analysis of bubbly and slug regimes 
of gas-liquid flows and is appropriate to the situation at 
hand. 


The limit of counter-current flow, called flooding as the 
pipe fills completely with liquid, predicted by the drift 
flux model can be seen on a plot of solutions of the 
momentum balance equations, (28) and (29), for constant 
liquid height in the drain pipe. Such a curve is plotted in 
figure 23, with superficial liquid velocities plotted 
horizontally and superficial vapor velocities plotted 
vertically, with each curve representing a different drain 
pipe liquid height. The region of negative liquid velocity 
and positive vapor velocity, called the limited region of 
counter current flow in reference Wl, is bounded by the axes 
and the dotted line known as the flooding line. The 
flooding line is the locus of points tangent to the j,,j, 
curves in the limited counter-current flow region, since an 
increase in the magnitude of either phase velocity beyond 
the flooding line is condition for which no steady flow 


Solution is possible and a flow transition must occur. 
Conclusions 


The suggested model, reflecting the influence of important 


physical features for processes influencing the measurement, 
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produces solutions that are stable and finite for both 
constant steam generator pressure and steam generator 
pressure gradients. Although a numerical error was made 
when solving the two phase flow equations (equations (28) 
and (29)), the steady state solutions appear to correspond 
well to the physical situation. The solutions produced for 
steam generator pressure ramps indicate that a potential 
shortcoming of the model is that some transients might cause 
the differential pressure between the condensing pot and the 
steam generator to exceed the values of differential 
pressure used to pre-compute superficial liquid and vapor 


velocities in the drain pipe. 
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Chapter 5 


Smart Instrument Applied to Steam 


Generator Level 


Having postulated a steam generator level measurement model 
capable of determining the pressure in the condensing pot 
and the liquid level in the drain pipe as functions of time, 
is appropriate to consider how successful a smart instrument 
including such a model might be at improving the accuracy of 
level measurement over a standard level sensing device. 
Since no operational plant level data is available for 
comparison, the Prism (reference Kl) program is used to 
generate steam plant data for an up power and a down power 


transient at high power. 


5.1 Sensor Dynamics, Information 


The relationship between steam generator water level and the 


differential pressure cell output is: 


MORD Ar ۲ ( OI) 0D) js Teu 


Tea (Tal -Py(Pag) 19 | Paol Tae) -Py (Pag) 


Equation (51) is the same as equation (2), except that the 
functional dependence of all the variables has been 
specified. Thus the liquid density, Pa. is a function of 
the temperature in the downcomer, the vapor density, p,, is 
a function of the steam generator pressure, and the 
reference leg liquid density, f,, is a function of the 


ambient temperature in the reactor building. 
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For a conventional level instrument, we can make some 
simplifications to make the level easier to compute. Over 
the pressure range usually encountered in commercial power 
reactors between 0 and 100 percent power, the vapor density 
changes from 36.5 kg/m? to 39.5 kg/m? so using an average 
value of 38.0 kg/m? for the vapor density introduces an 
error of no more than 4 percent. Assuming that the reactor 
building ambient temperature does not change, the reference 
leg density can be regarded as constant at 988 kg/m,. The 
conventional instrument can also neglect the pressure 


differential between the generator and the condensing pot. 


5.2 Downcomer Temperature Effects 


As noted in the description of steam generator internals, 
the steam generator downcomer is the region between tube 
sheet wrapper and the generator vessel, below the feedwater 
inlet ring. Prior to entering the tube bundle area, 
feedwater returning from the condenser mixes with saturated 
liquid that has been trapped by the moisture separators in 
the downcomer. The saturated liquid from the moisture 
separators also preheats the feedwater. The density of the 
downcomer liquid is related to its temperature. 
Unfortunately, downcomer temperature is not a measured plant 
variable, although feedwater temperature is. To estimate 


downcomer temperature, a heat balance is performed on the 


downcomer: 
Mac Cp Tae = Mey Cp Tet (Mic) Cpl sat (52) 
where hac â mass flow rate in rate in downcomer; 
Co 4 specific heat capacity of water; 
Tac 4 temperature of downcomer water; 
Me, â mass flow rate of feedwater; 
ee 4 temperature of feedwater; and 
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Tuc A saturation temperature of generator 
pressure. 
assuming that the specific heat capacities in equation (52) 
are approximately equal.  Rearranging (52): 


= Ihe T e, * une) a 


dc I (53) 


The feedwater and steam mass flow rates of the steam 
generator are measured plant variables, but the downcomer 
mass flow rate is not. The downcomer mass flow rate can be 
estimated using the Circulation Ratio. The Circulation 
Ration is a measure of the quantity of water entering the 
tube bundle section of the steam generator (downcomer flow) 
compared tot he total feedwater flow. Circulation ratio is 


defined as: 


Mass Flow Rate Entering the Tube Bundle _ Mac 


CR= 
Mass Flow Rate of Feedwater Mr 


wW 


Typically, the Circulation Ratio decreases from about 25 at 
lower powers to about 5 for high powers. By estimating the 
Circulation Ratio for the power level and multiplying it by 
the feedwater mass flow rate, it is possible to generate a 


reasonable value for downcomer mass flow rate. 


For the conventional instrument, it is assumed that a good 
estimate of downcomer temperature would be made and used to 
calculate a constant downcomer density. It turns out that 
the downcomer temperature is not a very sensitive function 
of feedwater temperature. For a feedwater flow rate of 20 
kg/s, feedwater temperature of 132°C, saturation temperature 
TT 290 °C, and a Circulation Ratio of 20, (corresponding to 


about 5 percent power), the downcomer temperature is 282°C. 
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For a feedwater flow rate of 400 kg/s, feedwater temperature 
5 209 مب‎ saturation temperature of 286°C, and a Circulation 
Ratio of 5 (corresponding to nearly 100 percent power), the 
downcomer temperature is 282?C. Clearly, the conventional 

instrument is capable of using a single value for downcomer 


liquid density, 747 kg/m? based on T,,-282?C. 


A smart level instrument could accept feedwater flow rate 
and feedwater temperature as inputs and, with knowledge of 
the details of Circulation Ratio as a function of power, 
calculate an instantaneous downcomer temperature using 
equation (53). This calculation is performed by the smart 


instrument in this report. 


Eo Inputs and Outputs 


Figures 24 and 25 show block diagrams for the conventional 
and smart instruments. The conventional instrument's only 
variable input would be the output from the differential 
pressure cell. The only output for the conventional 
instrument is level, L,. The smart instrument would have 
the differential pressure cell input as well as the 
temperature of the feedwater, T,,, the mass flow rate of the 


feedwater, ih,,, the Circulation Ratio appropriate for the 


power level, CR, the steam generator pressure, p,,, and the 
temperature of the reactor building, Tame. The smart 
instrument would provide two outputs: level, L, and a 
flooding warning to alert the operator of an approach to 


flooding in the drain pipe. 
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Equation (31) |— Level 





Figure 24: Conventional Level Instrument Block Diagram 


Equation (51) CE 





Figure 25: Smart Level Instrument Block Diagram 
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5.4 Calibration 


The output data from the Prism program after running 30 
minute simulations of up and down power transients ls used 
to calibrate the conventional and smart instruments. For 
the purposes of comparison, this procedure assumes that the 
steam generator data recorded by Prism represents the actual 
value or truth. The calibration procedure was to compute 
the downcomer water temperature as a function of time using 
equation (55), assuming a Circulation Ratio of 5 for the 
power level, to find the downcomer liquid density, Q, as a 
function of time, compute the vapor density, p,, as a 
function of time based on the Prism data for generator 
pressure, use fixed values of reference leg density, p,=988 
kg/m’, reference leg height, H-2 m, and use the Prism data 
for generator level to solve equation (1) for the 
differential pressure cell output, AP, as a function of 


eime. 


Naturally, using this calculation procedure, it is not 
instructive to compare the smart level instrument calculated 
level output with the level signal recorded by Prism. They 
are identical because the procedure and data used by the 
smart instrument to calculate steam generator level from 
differential pressure cell data is just the inverse of the 
procedure used to derive the differential pressure cell data 


from the Prism level information. 


SS Comparison 


The time varying differential pressure cell output, AP(t), 
extracted from Prism data by the procedure described in the 
previous section, was applied to the simplified calculation 
(51) performed by the conventional instrument and to the 


more detailed calculation performed by the smart instrument. 
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Comparison for Down Power 
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Figure 26: Comparison of Level Signals for Down Power 


Figure 26 is plot of the level signals as functions of time 
for the smart and conventional instruments for a down power 
from 100 to 90 percent power. Recall that the difference 
between the two calculated levels is that the conventional 
instrument uses an average value for vapor density, constant 
downcomer liquid density, and neglects the differential 
pressure between the condensing pot and steam generator. 
The smart instrument calculates the vapor density and liquid 
density at each time step, and does not neglect the pressure 
difference between steam generator and condensing pot. 
Figure 27 show that the level error between the conventional 
and smart instruments varies, but is generally below .5 
percent. 


85 





Level Error - Down Power 
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Figure 27: Level Error For Down Power 


Figure 28 is a comparison of the level signals for an up 
power transient from 90 to 100 percent. From figure 29, the 
minimum deviation between the level signals for the up power 
iS .4 percent, but at no time is the error greater than 1 


percent. 
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Figure 29: Level Error For Up Power 


Conclusions 


The operation and accuracy of a smart instrument was 
investigated. The model was used by the smart instrument 
was developed in chapter 4. The input variables for the 
smart instrument are steam generator pressure, feedwater 
flow rate and temperature, reactor building temperature, 
Circulation Ratio, and differential pressure cell output. 
For the basis of comparison, a "conventional" level 
instrument was developed by analyzing the behavior of the 
variables in equation (51) over the range of steam generator 


pressures reasonable during power operation and taking all 
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quantities in equation (51) constant except tor the 
differential pressure cell output. The only input variable 
to the conventional instrument is thus differential pressure 
cell output. Since not operational differential pressure 
cell data was available, the Prism program was used to 
generate level information as a function of time for two 
power transients and the differential pressure cell data for 
both transients was extracted from the Prism level data. To 
Ehe extent that the simulated differential pressure signal 
is valid, the smart instrument was more accurate than the 
conventional instrument, but the difference between the two 
level signals over the thirty minutes of simulation never 
exceeded one percent. It is likely that process noise, 
normal sensor fluctuations, and meter error would completely 


obscure such a small accuracy improvement. 


Although the smart instrument is not significantly more 
accurate than the postulated conventional level instrument 
over the range of pressure changes produced by the small 
power transients (ten percent) imposed on the Prism program, 
the smart level instrument is capable of estimating the 
liquid level in the drain pipe between condensing pot and 
Steam generator while the conventional instrument is not. 
Thus, the smart instrument, while being only slightly more 
accurate, 1S capable of identifying and warning the operator 
of a violation of the assumption common to both smart and 
conventional level instruments: that the reference leg 
height is constant. While a utility's response to a smart 
instrument reporting that drain pipe flooding during power 
operations would no doubt be exactly the same as if a 
conventional level instrument's output produced unexpected 
jumps in level characteristic of drain pipe flooding (i.e., 
replace the drain pipe or condensing pot with different size 
piping), the smart instrument's ability to infer variables 


not easily measured and take action based on that 
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tion has been demonstrated. 
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Chapter 6 


Conclusions 


A procedure for designing and implementing a smart 
instrument was applied to a Pressurized Water Reactor steam 
generator water level measurement. Important physical 
effects on generator level were considered and applied to 
find an overall model capable of inferring desired variables 
that are not directly measured. The model is stable for 
constant steam generator pressure and is also stable for a 
limited range of steam generator pressure transients. A 
conventional level instrument, using a simplified procedure 
to calculate level, was chosen for comparison with the smart 
level instrument. Both instruments received differential 
pressure cell data extracted from Prism level for two 
transients. As expected for this comparison, the smart 
instrument output conformed more closely to the "true" steam 
generator level, but the difference between the smart 
instrument and the conventional instrument was never more 


than 1 percent. 


۲ Significance of Corrections 


To understand why there was so little difference between the 
smart instrument level data and the conventional instrument 
data, it is necessary to compare the magnitude of the 
corrections made by the smart instrument to the differential 


pressure signal. 


The differential pressure cell output extracted from Prism 


output data is on the order of 365 kPa. The maximum steam 


I 





| 


| generator to condensing pot pressure difference calculated 


was 1 kPa. Clearly, it is safe to conclude that the 


differential pressure is negligible. 


The smart instrument calculates the downcomer temperature 
and uses this temperature to adjust the liquid density in 
the downcomer at each time step. The conventional 
instrument uses a constant downcomer temperature of 282°C. 
The range of downcomer temperatures calculated by the smart 
instrument is 280?C (0,2750 kg/m?) to 282°C (p,=747 kg/m’). 
This small variation in downcomer temperature produces a 


variation in liquid density of only 0.49 percent. 


The steam generator pressure changes for both the up and 


down powers were on the order of 100 kPa from a nominal 


Steam generator pressure of 7 MPa. This pressure deviation 


results in a maximum vapor density difference over the 


entire maneuver of 1.6 percent. 


Originally, the temperature of the reactor building was 
intended to be an input to the smart instrument since it 
seemed plausible that a steam leak in the reactor building 
would significantly raise the ambient temperature in the 
vicinity of the reference leg, thus decreasing the density 
of the liquid in the reference leg. Indeed, a steam leak 
that raised the reactor building temperature from 50°C to 
200°C would change the density of the reference leg from 988 
kg/m? to 865 kg/m’, a change of 12.5 percent. However, the 
magnitude of the differential pressure cell output, AP in 
equation (51), dominates the value of computed level so the 
resulting effect on generator level of a reactor building 
temperature change of this magnitude would be only .7 
percent. A steam leak into the reactor building would 
significantly raise the ambient humidity, making all 


instrumentation suspect of erratic operation and place the 
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plant at risk of uncovering the core. Under the 
circumstances, a .7 percent change in indicated steam 
generator level would not be as important as the other 


effects on the reactor plant. 


Obviously, none of the simplifications and assumptions made 
by the conventional instrument envisioned for comparison 
with the smart instrument produce an error in excess of 1.6 
percent. This small deviation from the "true" steam 
generator level is really too small for the operators to 
notice, so the smart instrument measuring steam generator 
level does not provide a realistic increase in indication 
accuracy. It still would be useful for warning the operator 


of the approach to a flooding situation in the drain pipe. 


-6.2 Smart Instrument Design 


This research indicates there are several considerations to 


observe when designing a smart instrument: 
-select the variables for measurement; 


«choose a model that incorporates as much detail about 
the physical process as necessary to achieve the 
measurement accuracy desired and the computational 


resources of available microprocessors; 
«evaluate the model stability and accuracy of results; 


«Compare the smart instrument with a conventional 
instrument to evaluate the magnitude of the accuracy 


improvement; 


«compare information available to the operator using 


the conventional instrument to that provided by the 
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smart instrument, even if not significantly more 
accurate, the designer may still opt for the smart 
instrument because of the extra information it 


provided. 


6.3 Flooding 


The suggested model describing steam generator level 
dynamics, to the extent it closely approximates the actual 
variables involved, is capable of calculating the liquid 
level in the drain pipe between the steam generator and the 
condensing pot directly. By observing when the liquid level 
reaches or almost reaches the top of the drain pipe, the 
smart level instrument envisioned by this research could 
warn the operator that its level output data is on the verge 
of unreliability. It is important to note, that the ability 
of the smart instrument to compute a solution of equation 
(49) for the state variables is predicated on the condition 
that flooding in the drain pipe does not occur, since 
equations (28) and (29), used to relate the liquid height in 
the drain pipe and the pressure gradient between the 
condensing pot and steam generator to the superficial 
velocities in the pipe, are not valid after flooding occurs. 
The pertinent question then is how does the smart instrument 
recover from just such a flooding condition since its model 


is now no longer valid? 


To recover from a flooding condition, the smart instrument 
should be programmed with a routine to estimate the steam 
generator water level based on the reduced reference leg 
height discussed in section 4.1 and depicted in figure 5. 
The instrument should warn the operator that flooding has 
occurred, but still be able to present its best estimate for 
water level. Next, the routine used to find the superficial 


velocities from the tables that have been pre-computed must 
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be capable of responding to the case when the input 
arguments, liquid height in the drain pipe and pressure 
gradient, fall outside the bounds of the tabular data. The 
smart instrument should then try to calculate level after 
discarding its prior data, in effect re-zeroing itself. It 
should continue to zero itself and try to calculate 
transient values of the state variables until the a the 
drain pipe liquid level is reduced and the pipe is no longer 


flooded. 
Conclusions 


This research discusses how and why smart instruments might 
be applied to nuclear power plants. The features that 
distinguish a smart instrument from a conventional 
instrument are a process model and the ability of the smart 
instrument to make decisions based on plant variables and 
model output. A good model of the physical process enables 
the instrument to monitor measurement results and warn the 
user if the data received is not consistent with redundant 
sensors or trends indicated by sensors for different 
variables. The important issues to consider when designing 
a smart instrument were explored and actually illustrated by 
a specific example: design of a smart instrument for steam 
generator water level measurement. A model of steam 
generator level variations as a function of power was 
proposed to reduce the number of assumptions necessary for 
level computation and to infer unmeasured variables germane 
to the reliability of the remaining assumptions. The Prism 
program was used to generate simulated level data, from 
which an estimated differential pressure cell output signal 
was extracted. The differential pressure cell signal and 
other Prism data were input to the smart level instrument. 
The "smart" level produced was compared to a conventional 


level instrument’s output. The smart level signal was only 
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marginally more accurate than the conventional level signal 
but the smart instrument is capable of inferring important 
variables that the conventional instrument cannot. Finally, 
the smart instrument must be capable of recovering from 
brief operation outside the limits of its model, otherwise 
the designer risks that certain transients could render the 


instrument's subsequent output unreliable. 
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Chapter 7 


Recommendations for Future Work 


There has not been sufficient time during the term of this 
research to pursue all the issues involved in smart 
instrumentation implementation as thoroughly as desired. 
Topics related to a smart level instrument worthy of 
additional investigation include: determining the criteria 
affecting the theoretical stability of the one-step, 
implicit difference equation (equation (48)) used to compute 
the solution to the linearized equations, investigating 
model behavior at lower steam generator pressures, 
numerically analyzing of the solution method, and obtaining 
experimental data on flow transitions and pipe liquid height 
for two phase counter-current flow and comparing the 
experimental results to the predicted results. Applying the 
design methodology suggested in this research or another 
method to the measurement and display of other nuclear 
reactor plant variables, such as reactivity measurement and 
control, would be of interest for the continued study of 


Smart instruments in general. 


7.1 Theoretical Stability 


Before implementing the proposed model for steam generator 
level measurement, the specific conditions for stability of 
equation (24), the general differential equation, and 
equation (48) the one-step implicit difference, should be 


determined. 


7.2 Operation at Lower Pressures 
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The steam generator level model developed in this research 
was only applied to level measurement at normal steam 
generator operating pressure, about 7.0 to 7.6 MPa. For the 
smart instrument to really replace conventional level 
instruments, it must be capable of operating over the entire 
range of generator pressures encountered in practice, 
extending down to atmospheric pressure. To investigate 
model performance at other pressures, solutions to the 
momentum balance equations, (28) and (29), must be computed 
for several pressures between 7.6 MPa and atmospheric since 
the liquid and vapor densities, taken as constants when 
iteratively solving equations (28) and (29), are strong 


functions of pressure in the steam generator. 


7.3 Physical Testing 


As the comparison of results in Chapter 5 illustrated, the 
smart level instrument employing a detailed model was only 
marginally more accurate than a conventional level 
instrument using a simplified model. Still, the smart 
instrument is capable of calculating the transient liquid 
level in the drain pipe and can therefore detect the onset 
of flooding in the drain pipe and subsequent erroneous level 
indication. Since the ability to predict the drain pipe 
liquid level is the distinguishing feature of the smart 
instrument, the accuracy of the prediction should be 
evaluated by physical experiment. The physical data 


obtained could be used to improve the model used as well. 


7.4 Numerical Analysis 


Computing the solution to equation (48) is very time 
consuming. The solution for the positive pressure ramps 
requires several hours to complete on a DECstation 3100, a 


particularly fast workstation. It would be potentially 
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beneficial to further analyze the equations involved and 
construct an algorithm designed to enhance the speed of the 
solution. Certainly, the speed of the calculations could be 
significantly enhanced by the use of a compiled high level 


computer language for the computations. 


7.5 Application to Other Variables 


The application of smart instruments to measurement and 
control of other nuclear power plant variables should be 
investigated. In addition to evaluating and displaying the 
value of a simple plant variable, use with complicated 
Systems like reactivity or steam generator water level 
control should be explored. A smart instrument could be a 
key component of an expert system used for power plant 
operation and casualty control. The smart instrument could 
be responsible for redundancy and consistency checks while 
an external computer evaluates the over all system, 
postulates probable outcomes to certain actions, and can 
recommend corrective action. Specific applications 
utilizing the ability of smart instruments to communicate 
with each other should also be explored. Sensors could be 
networked together and either polled by a control system or 


signal the control system when they have data to transmit. 
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Appendix 


The PRO-MATLAB program, numerical processing software based 
on LINPACK and EISPACK codes, is used to solve equation (48) 
for each of the postulated transients.  PRO-MATLAB is an 
optimized program that is well-suited for high speed 
calculations on large matrices and vectors. The program is 
capable of processing commands that have been saved to a 
text file, called a macro, and allows the user to create his 
own functions. Neither functions or macros are compiled, 
but they were used to develop prototype algorithms to 
compute solutions to equation (48). The same routines could 
be compiled and used by the smart level instrument to 
perform similar calculations. The matlab macros and 


functions used are included in the following pages. 
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$macro to start calculating a solution to equation (48) 
 stmdat 

gets rhol list, rhov list, Hl list, Hv list, pressure 
variables 5 

gets table h and table press values (for jv and jl data), 
jvdat and jldat 

dtz1*.5; 

L=2.5;%length of drain pipe 

Dt=.0254;%drain pipe dia 

Dc=2* .0254; $condensing pot dia 

Me=3.5*2.54e-2;%length of condensing pot 

dpdx=-73 .27493046894671; 

startpsg=7.34e6; 

endpsg=7 .0400e6; % 

etoptrans=225/dt; 

minutes-4;£number of minutes to simulate 
Enput-[startpsgtdpdx*L;.71227377;0.00546943]; 


deltas-[-30;.001;.0005];$deltas to use in computing 
derivatives 

E =[]; 

ul]; 

prsg=[]; 

$X is the set of intermediate values of the result for each 
time $step 

for n=-5/dt:minutes*60/dt 

%gcontrols psg (t)-> 

if n<=0 

psg=startpsg;%initial value 

elseif n<stoptrans+1 

psg=startpsg-n* (startpsg-endpsg) /stoptrans;%for ea time 
$step during power chg 

else 

psg=endpsg; 

end 
X -[X input]; 
prsg-[prsg; psg]; 
psg 
% display dp/dx and dp/dx + dp 

(input (1) -psg) /L 

(input (1) —-psg) /L+deltas (1) 

[X, f£]=solveit (psg, 70000, dt, input, L,Dc,Dt, Lc, deltas, ... 
EE c wyidat,rhol list,rhov list,Hl list,Hv list,.. 
Mabile press,table h, pressure); 

input=X; 

disp('finished time step') 

n*dt 

EI f]; 
end 
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$this is the macro stmdat 

$get steam table data in flat file form 

load stmdata.prn 

n-16;$number of pressure data pairs 
pressure-stmdata(1:n); 

Boel list=ones(1l:n)’ ./stmdata(nt1:2*n) ; 

Env Irst-ones(1:n)'./stmdata (2*n*1:3*n); 

Emu ust-stmdata(3*nt*1:4*n); 

fee ttst—stmdata (4*n+1:5*n) ; 

%load outdata3.prn 

load outdata4.prn 

nl-41;$number of pressure data pairs 

teen ۱0۶ ۰0012 ۰002152 .002; .0025; .003; .0035; .004; .008;... 
۱۳ 016; .020; .023;025];%order of liquid heights 


table press=[{600;570;540;510; 480; 450; 420; 390; 360; 330;300;... 
0 310—210;180;150;120;90;60;30;0;-30;-60;-90.1;-120.1;... 
meee —7180.2;-210.2;-240.2;-270.2; -300.3; -330.3;-360.3;... 


iE —420.4;—450.4;—480.4; -510.4; 540.5; —570.5;-600.5]; 
jvdatzoutdata4(1:14*n1); 


jldat=outdataf4 (14*n1+1:28*n1); 
end 
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Rumnetion [X,f] moO LVelE( psa, dosgq,dt, inputs, by, Oc, Pe, uc, ... 
ENUvdat,3ldat,rhol ist roy list; Hise, av ast, 

table press,table h,pressure) ; 

Sfunction to solve the equation and compile the answers 


$psg is initial generator pressure 

$dpsg is delta generator pressure 

$dt is time interval of manuever 

sinputx is inputx-[pc; alpha erh] 

$L is length of drain pipe 

$Dc is condensing pot diameter 

$Dt is drain pipe diameter 

$a is deltas of Xx variables for numerical differentiation 
۳۳۳ ۱1۳۴۶ 1۵1۳۴12 c;dh] [-30;.0005;.0005] 

Eusget jvdat,jldat, rhol list, rhov list, Hl list, 

Hv list,table press, 

$table . h from running stmdat 

$pressure is the list of pressures for which steam table 
data is available in 

$the file stmdata.prn read by stmdat.m 


$define given deltas to use for numerical derivatives 
dp=d (1) ; 

da=d (2) ; 

dh=d (3); 

5 

$for n=1:5 


$ psg=psg-n*dpsg/600;%for each time step 
$define state variables 

alpha c-inputx (2); 

pc=inputx (1); 

hl=inputx (3); 

$now build the df/dx vector 
Ess—rrrst(alpha c,psg,pc,hl,L,Doc,Dt,dp,Jjvdat,;jldat, 
table press,table h); 
Bee 2—Second(alpha c,psg,pc,hl,L,Dc,Dt,da, jvdat, jldat, 
table press,table h); 
Esr-thrrd(alpha c,psg,pc,hl,L,Dc,Dt,dh,jvdat, jldat, 
Fable press,table h); 
fee fourth (alpha c,psg,pc,hl,L,Dc,Dt,dp, jvdat, jldat, 
pable press,table nh); 
Eus2—rrfth(alpha c,psg,pc,hl,L,Dc,Dt,da,jvdat,jldat, 
Fable press,table EST 
af23= =sixth (alpha c, Ps PENI; L De, Dt dn i vdat; Jidat, 
table press,table h); 
df31=seventh (alpha_c,psg,pc,hl,L,bc,Dt,dp, jvdat, jldat, 
Fable press,table _h); 
Hso2-erghth(alpha c,psg,pc,hl,L,Dc,Dt,da, jvdat, jldat, 
Fable press,table h); 
df33=ninth (alpha_c,psg,pc,hl,L,Dc,Dt,dh, jvdat, jldat, 
Fable press, table h); 
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ren dfl2 dfl3;df21 df22 df23;df31 df32 df33]; 


$ 


get the f vector 

f-findf(alpha c,psg,pc,hl,Dc,Dt,L,jvdat,Jjldat,... 
table press,table Rh) 

% 


$get the A vector 
% 


ELUOMUIOA (pc,alpha c,hl,rhol list,rhov list,Hl list,... 


Hv l1ist,Dc,Dt,L,Lc,pressure); 
X=inv((A/dt - df))*f + inputx; 
disp(’ finished time step’) 


end 
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emnetion [dfldpc]=first 
(alpha c,psg,pc,hl,L,Dc,Dt,dp,jvdat,jldat,press,givenh) 
% 


function to compute the derivative of fl wrt pc 

% 

$functions first.m through ninth.m use similar arguments: 

$ 

$alpha c is the vapor volume fraction in the condensing pot 
$psg is the steam generator pressure 

spc is the condensing pot pressure 

$hl is the liquid height in the drain pipe 

$L is the length of the drain pipe 

Dc is the diameter of the condensing pot (a horiz cylinder) 
$Dt is the diameter of the drain pipe 

$dp or da or dh are the deltas used to compute the 
derivatives Swrt 

Ecspha c, or hl 

$press and givenh are the P’ and h data used to compute the 
jvdat %and jldat 

$tables 


theta=.0164; 
% 


$get steam data 
ExbsTsat,mu,k,rhovl,rhol] 
Eb UTSsat,mu,k,rhov2,rhol] 
% 


SEMIRNEO (PCJ; 
SEMINEO (DEFOD); 


$find mdot 

% 

mdotwl = findmdot (alpha _c,Dc,Dt,pc,theta); 

mdotw2 = mdotwl;%same as 

mdotw2=findmdot (alpha_c,Dc,Dt,pc+tdp,theta); 

$ 

nd gv 

$ 

= تا 

findjvjl(hl, (pc-psg) /L, givenh, press, jvdat, jldat); 
[j12 jv2] = 

findjvj1(hl, (pc-psg) /L+dp, givenh, press, jvdat, jldat) ; 
$ 


A-(Dt/2)^2*pi;$area of pipe 

$fl-rhovl1*jvl1*A-mdotw1 

dfldpe = ((rhov2*jv2*A - mdotw2)-(rhovl*jvl1*A - mdotw1))/dp; 
end 
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۳۳۲۰۲۰ 1 5۳7 ۱04۴ 1026[ 2566010 )81۳۳۵ c,psg,pc,hl,L,Dc,... 
Dt,da, jvdat, jldat,press,givenh) 
$ 


Seunction to compute the derivative of fl wrt alpha _c 
$ 

¿function [dfldac]) = second 

Espha c,psg,pc,hl,Dc,Dt,da,jvdat, jldat) 

$ 

theta=.0164; 


$get steam data 


[a,b,Tsat,mu,k,rhovl,rhol] = stminfo (pc); 

% 

$find mdot 

$ 

mdotwl - findmdot (alpha c,Dc,Dt,pc,theta); 
mdotw2 - findmdot (alpha ctda,Dc,Dt,pc,theta); 
$ 

$find jv 

$ 

EI v1) = 

findjvjl (hl, (pc-psg) /L,givenh, press, jvdat, jldat); 
۲۱۱2 v2] = 


Enndjv3l (hl, (pc-psg) /L,givenh,press,Jjvdat, jldat); 
$ 


A-(Dt/2)^2*pi;$area of pipe 

$f1=rhov1*3v1*A-mdotwl; 

dfldac - ((rhovl*jv2*A - mdotw2)-(rhovl*jv1*A - mdotw1))/da; 
end 


LIU 


EU ceron[dfldh]-third(alpha c,psg,pc,hl,L,Dc,... 
Boah, jvdat, jldat,press,givenh) 
% 


$function to compute the derivative of fl wrt h 
2 


o 


Eencetion [dfldpc] = third 
(alpha c,psg,pc,hl,L,Dc,Dt,dh, jvdat, jldat) 


theta=.0164; 


% 

%get steam data 

[a,b, Tsat, mu, k, rhov1, rhol] = stminfo (pc); 
% 

$find mdot 

$ 

11001-171 = findmdot (alpha _c,Dc,Dt,pc,theta); 
$ 

5 ۶ 1110 jv 

% 

۱ 1 10۳11 = 

Endjv3l (hl, (pc-psg) /L,givenh,press,Jjvdat,jldat); 
[j12 jv2] - 


ENUudgjvjl (hltdh, (pc-psg) /L, givenh, press, jvdat, jldat); 
$ 
A-(Dt/2)^2*pi;$area of pipe 


dfldh 2 ((rhovl*jv2*A - mdotwl)-(rhovl*jv1*A - mdotw1))/dh; 
end 


Ti 


ED :vonidf2dpc]-fourth(alpha c,psg,pc,hl,L,Dc,Dt,... 
dp, jvdat, jldat, press, givenh) 
% 


$function to compute the derivative of f2 wrt pc 
% 

Eunetron [dfldpc] » fourth 

Emo u33sc.psg,pc,hl,Dc,Dt,dp, jvdat, jldat) 


theta=.0164; 
% 


$get steam data 

MAY], Tsat,mu,k,rhovl,rhol] 
NAO, Tsat,mu,k, rhov2,rhol] 
% 

5 111 1 

$ 

mdotwl = findmdot (alpha_c,Dc,Dt,pc,theta) ; 

mdotw2 = mdotwl;% same as 
mdotw2=findmdot (alpha _c,Dc,Dt,pctdp, theta) ; 

% 

$find jv 

% 

11 jvi] = 

findjvjl (hl, (pc-psg) /L, givenh, press, jvdat, jldat); 

] 12 jv2] - 

findjvjl (hl, (pc-psg) /L+dp, givenh, press, jvdat, jldat) ; 
% 


stminfo (pc); 
stmamnito (peddpb); 


$findQ 

% 

Ql=findQ (pc) ; 

Q2=findQ (pctdp) ; 

A-(Dt/2)^2*pi;$area of pipe 
f2-rhovl*jvl*Hvl*A-mdotwl*H11-01; 

df2dpc - 

NE Ucu2593v2*Hv2*A-mdotw2*H12-02) -(rhovl*Jjvl*Wvl*A-mdotwl*Hll 
01))/dp; 

end 


1I2 


Nux cuonun[idf£2dac]-fifth(alpha c,psg,pc,hl,L,Dc,Dt,... 
da, jvdat, jldat, press, givenh) 
% 


Beeretioen to compute the derivative of f2 wrt alpha c 
% 


theta=.0164; 
% 


%get steam data 

NAL, Tsat,mu,k, rhov1,rhol] = stminfo (pc); 
% 

%find mdot 

$ 


mdotwl - findmdot (alpha c,Dc,Dt,pc,theta); 


mdotw2 Einomdst (alpha ctda,Dc,Dt,pc,theta); 

% 

%find jv 

% 

lI vl] = 

findjvjl(hl, (pc-psg) /L,givenh,press, jvdat, jldat); 
% 

۶ ۶ 1 (0 


Ql-findQ (pc); 


A-(Dt/2)^2*pi;$area of pipe 
$f£2=rhov1*3v1*Hv1*A-mdotw1*H11-Q1 

df2dac - 

((rhovl*Jjvl1*Hvl*A-mdotw2*Hll-Q1)- (rhovl*jvl*Hvl*A-mdotwl*Hll 
Q1))/da; 

end 
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EUcconfífdf2dh]2sixthí(alpha C,psg,pc,hl,L,Dc,Dt,... 
dh, jvdat, jldat,press,givenh) 
% 


%function to compute the derivative of f2 wrt h 
% 

$function fat lan] s third 

Ey o»najc,psg,pc,hl,Dc,Dt,dh,Jjvdat, jldat) 


theta=.0164; 
$ 


$get steam data 

Ml, Tsat,mu,k,rhovl,rhol] = stminfo (pc); 
% 

%find mdot 

% 

mdotwl = findmdot (alpha_c,Dc,Dt,pc,theta); 


$ 
$find jv 
9; 


[jll jvi] = 

findjvjl (hl, (pe-psg) /L, givenh, press, jvdat, jldat) ; 
[j12 jv2] = 

findjvjl (hl+dh, (pc-psg) /L, givenh, press, jvdat, jldat) ; 
% 


Efrindo 

% 

Ql=findQ (pc) ; 

A=(Dt/2)*2*pi;%area of pipe 
f2=rhovl*jv1l*Hv1l*A-mdotwl1*H11-Q1; 

df2dh = 
((rhovl*jv2*Hvl*A-mdotwl1*H11-Q1)-(rhovl1*jvl1*Hvl*A-mdotwl*Hl1 


-91) ) /dh; 
end 
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£unction[df3dpc]=seventh (alpha c,psg,pc,hl1,L,Dc,Dt,... 
dp, jvdat, jldat,press,givenh) 

% 

$function to compute the derivative of £3 wrt pc 

% 


theta=.0164; 
% 


$get steam data 

۱۶ bTSat,mu,k,rhovi,rhol11] 
Bon, Tsat,mu,k,rhov2,rhol2] 
% 


stminto (Pe); 
stminfo (pc+tdp); 


$find mdot 

$ 

maorwl = findmdot (alpha c,Dc,Dt,pc,theta) ; 

mdotw2 = findmdot (alpha _c,Dc,Dt,pctdp, theta) ; 

$ 

$find jv 

$ 

[j11 jv1] = 

findjvjl (hl, (pe-psg) /L, givenh, press, jvdat, jldat) ; 
[312 jv2] - 

Sanajvjl (hl, (pc-psg) /L+dp, givenh, press, jvdat, jldat) ; 
$ 


A-(Dt/2)^2*pi;$area of pipe 

f3=rholl*3jJ11*A+tmdotwl; 

df3dpe = ((rhol2*j1l2*A + mdotw2)-(rholl*jli*A + mdotwl))/dp; 
end 
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jvdat, jldat,press,givenh) 
% 


Seanetion to compute the derivative of f3 wrt alpha c 
+ 


theta=.0164; 
$ 


$get steam data 

ExbUrsat,mu,k,rhovi,rholl] = stminfo (pc); 
$ 

$find mdot 

$ 
mdotwl 
mdotw2 
$ 
$find jv 

$ 

[j11 jv1] = 

Sima jvjl (hl, (pc-psg)/L,givenh,press,]jvdat,jldat); 
Eu 7jv2] = 

Eindsv3i(hl, (pc-psg) /L, givenh, press, jvdat, jldat) ; 
$ 


amador (alpha c,Dc,Dt,pe, theta); 
findmdot (alpha c+da,Dc,Dt,pc,theta) ; 


A-(Dt/2)^2*pi;$area of pipe 

$f3=rho11*311*A+mdotw1; 

df3dac = ((rholl*j12*A + mdotw2)-(rholl*jl1*A + mdotwl))/da; 
end 
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Ncrson[idf3dn]enznth (alpha c,psg,pc,hl,L,Dc,Dt,... 
Envdat,jldat,press,givenh) 
$ 


$function to compute the derivative of f3 wrt h 
$ 


theta=.0164; 
% 


5061۲ steam data 

IND Tsat,mu,k,rhovl,rholl] - stminfo (pc); 

$ 

%gfind mdot 

% 

mdotwl = findmdot (alpha_c,Dc,Dt,pc,theta); 

$ 

$find jv 

$ 

[j11 jv1] = 

Minds) (bl, (pc-psg) /L,givenh,press, jvdat, jldat) ; 
[312 Jv2] = 

findjvjl (hl+dh, (pc-psg) /L, givenh, press, jvdat, jldat) ; 
% 

A=(Dt/2)*2*pi;%area of pipe 


df3dh - ((rholl*jl2*A - mdotwl)- (rholl*jll*A -* mdotw1))/dh; 
end 


2:157 





۰ ۲۶۰۰ 1ظ ,معط ,قوجر م6 21۳8۳۵) 221۳0۶[ ۱۶ 11 ۲ ۳9 
L, jvdat, jldat,press,givenh)‏ 

$function to compute the f vector 

$ 

Eha € is the vapor volume fraction in the condensing pot 
%psg is steam generator pressure 

$pc is condensing pot pressure 

$hl is liquid height in drain pipe 

$Dc is condensing pot diameter (a horiz cylinder) 

$Dt is drain line diameter 

$L is length of drain line 

$press and givenh are the P' and h data used to compute the 
jvdat $and jldat 

%tables 


$constants 

theta=.0146; 

$ 

$get steam data 

[Hl,Hv,Tsat,mu,k,rho_v,rho 1]=stminfo (pc); 
% 


$find mdot 


+ 


mdotw - findmdot (alpha c,Dc,Dt,pc,theta); 

$ 

۱ 11 jv, j1 

Pol gv) = findjvjl (hl, (pe-psg) /L, givenh,press, jvdat, jldat); 
$ 


%find Q 
Q-findQ (pc); 
$ 


A-(Dt/2)^2*pi;$area of pipe 
% 

%find f1 

fl-rho v*jv*A-mdotw; 

$ 

۳۳ 15:0 £2 

f2-rho v*jv*Hv*A-mdotw*Hl-Q; 
$ 

$find f3 
£3=rho_1*31*A+mdotw; 

$ 

> ۱۶12 ۶22 3 [ ۶ 

end 
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function[A]- buildA(pc,alpha c,hl,rhol ۱ e 
ist, Hv list,Dc,Dt,Lt,Lc,pssr) = 

% 

$function that builds the A matrix, xxx list are the stm 
table %lists of 5 

$data for the pressures read in from stmdata.prn by stmdat.m 
$pssr is the range of pressures used to extract steam table 
data 

$ 

$pc is pressure in condensing pot 

$alpha c is vapor volume fraction in condensing pot 

$hl is liquid hight in the drain pipe 

$Lt is the length of the drain pipe 

$Lc is the length of the condensing pot (assumed to be a 
horiz %cylinder) 

+ 

$uses: 1. stminfo.m 

% 2. deriv.m 

% 


Vc-(Dc/2)^2*pi*Lc;$volume of condensing pot 

[Hl Hv Tsat mu k rhov rhol]-2stminfo(pc);$stm table info 
$ 

$note: deriv(pssr,rhol list,pc,.01) computes the partial 
deriv of %rhol wrt pc at pc 


$compute all 

Bee zalpha c) *deriv(pssr,rhol list,pc, .Ol)+talpha_c*... 
(deriv(pssr,rhov_list,pc,.0Ol)))*Ve; 

% 

$compute a12 

al2=(rhov-rhol) *Vc; 

% 

$compute a21 

pr arpha c)*deriv(pssr,rhol list.*EHl list,pc,.01)-*... 
Enc deriv(pssr,rhov list.*Hv list,pc,.01))-1)*Vc; 

% 


$compute a22 

a22=(rhov*Hv -—- rhol*Hl1)*Vc; 

$ 

$compute a33 
a33=Dt/2*Lt*rhol* (sin (acos (1-2*h1/Dt)))72/sqrt (h1/Dt* (1-h1/D 
E); 

$ 

complete A 

A-[all,al2,0;a21,a22,0;0,0,a33]; 

end 


SIE 


function y - deriv(x,y,input,delta) 

% 

$computes the derivative of y wrt x at the input value OLI, 
using an x step 

$of delta 

$format is y = deriv(x,y,input,delta) 

y = (spline (x,y,input+tdelta)-spline (x,y,input))/delta; 


end 
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function [hl, hv, Tsat,mu,k,rho_v,rho 1] = stminfo (p) 

۱ ton to compute the properties of water given pressure 
E.. pressure in the range 1.5MPa to 9Mpa 

E... [hl,hv,Tsat,mu,k,rho v,rho 1] = stminfo (p) 


k-626e—-6;$thermal conductivity 

$steam table information 

press = [1.5e6;2e6;2.5e6;3e6;3.5e6;4e6;4.5e6;5e6;5.5e6;... 
6e6;6.5e6;7e6;7.5e6;8e6;8.5e6;09e6]; 


۳ ۶۱ 999. ۵0262۶ 549. 90 65 : 835 . 421 9: ٩22 .1 656۶7 10.044 ۰ 
8M 5557 787.9600; 777.7259; 767.8722; 758.3226; 749.0637; ... 
NA A 731.1545; 722.4390;713.7759; 705.2684] ; 


m X 11 -5955;10.0462; 12.6855; 15.0083; 17.5346;... 
E 722.7066; 25.3614; 28.0662; 30.8261; 33.6485;... 
E5525; 39.4836; 42.5080; 45.4339; 48.7924]; 


H1=1.0e+03*[.8447;0.9086;0.9618;1.0084;1.0498;1.0874;... 
8M 1545;1.1849;1.2137;1.2417;1.2674;1.2927;1.3171;... 
E7043 1.3637]; 
BE uUETFO3*[2.7899;2.7972; 2.8009; 2.8023;2.8020;2.8003;... 
8E 027 7942;2.7899;2.7850;2.7795;2.77735; 2.77669; 2.7599; ... 
2.7523;2.7446]; 
TD Sg. 29;212.37;223.91;233.84;242.54;250.33;257.41;... 
۱۳۳۱۱ ۱۱۱209 م9‎ 27/5. 55 280.827 285. 797 290.57 294.977 299.14۶ 303.3 
1]; 
$interpolate to get precise T and P in generator for given 
power %level 
$i and itl are the indices in the list that bracket the 
value %sought 
i-comp (p,press); 
if i--length(press) 

rho v-rhov (i); 

rho l-rhol (1i); 

hvzHv (i); 

hl=H1 (i); 

Tsat=TSAT (1); 
else 

۰۳ ۷ (p-press (1))./ (press (1+1)-press (1)) .* (rhov(i+1)-... 
rhov (1))+rhov (1) ; 

Ele I=(p press (i) ) ./ (press (i1+1) press (i) ) .* (rhol (1+1)~..- 
rhol (i) ) +rhol (i) ; 


hv-(p-press(i))./(press(itl)-press(i)).*(Hv(it1)-Hv(i))*... 
B (1); 


hl-(p-press(i))./(press(itl)-press(i)).*(Hl(it1)-Hl(i))t... 
B1(1); 
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We pepress(i))./ (press (itl)-press (i)) .*(TSAT(i+t1)-... 
ESAT(1))-TTSAT (i1); 
end 
mu-visc(Tsat,rho 1); 
end 
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function [index]=comp (value, list) 

%function index=compare (value, list) 

$function that compares each element of value with all the 
values in a 

$sorted list and returns the a vector containing the index 
of the closest 

%number below each element of value in list 


index-zeros(1,length(value)); 
for j=l:length (value) ; 
fom t—1l:length(list); 
if value (Jj) >=list (i) 
index (j) =1; 


end 
end 
if index (3) ==0 
disp(’’) 
disp('number is not in given range’) 
end 
end 
end 
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E meltion [0] =findQ (Pp) 

S$function to compute heat transfer rate as a function of 
pssr in $condensing pot 

S$function [Q]-2findQ (p) 


%constants 

g=9.807; 

Tamb=50; 

$wtd length. length is 3.5 inches half of which is 1.75 
inches=5.08e-2 

123.5/2*2.54e-2; 

$1=3.8le-2; 

%rin=7.42e-2; 

Srout=7.62e-2; 

rin=1.939*2.54e-2/2;%based on sched 80, 2 in diameter 
rout=2.157*2.54e-2/2; 

ksteel=50; 

hout=11; 

% 

$get steam data 
[hl,hv,Tsat,mu,k,rho v,rho 1]sstminfo(p); 

% 

<hfg=hv-hl; 

$compute hbar 

IN ie rho 1.^2*g*1^3./ (4A*mu*k.* (Tsat-Tsurf))).^.25; 
$hbar-24/3*k*1e3/l1*Nu;$1e3 needed since k is in kW 
‘compute UA 

SUA=((1/ (2*pix*rin*l*hbar) )+(log(rout/rin) / (2*pixl*ksteel))+. 


BEN ZX*pi*rout*l*... 

pbout)))^(-1); 

UA=((log(rout/rin)/ (2*pi*l*ksteel))+(1/ (2*pi*rout*1*... 
Bout)))^(-1); 

$compute Q 

Q-UA.*(Tsat-Tamb);$mass flow rate 

end 
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function [mdotw] - findmdot (alpha c,Dc,Dt,pc,theta) 
2 ze 

$function to compute the mass flow rate in the weir 
$ 

$alpha c is the vapor volume fraction in the condensing pot 
$Dc is the diameter of the condensing pot 

$Dt is the diameter of the drain pipe 

$pc is the pressure in the condensing pot 

$theta is the elevation angle of the drain pipe 

$ 

$uses: 

$ l.stminfo.m 

$ 2.interp.m 


$constants 

g=9.807; 

al=2; 

$find the liquid density 

[Ebr d.e,f,rhol] = stminfo (pc); 

$find the phi c that corresponds to the given alpha c 


Ec-nterp(alpha c);$call special routine to find phi c 


hlc-Dc/2*(1-cos(phi c));$1iq ht in condensate pot 
$ 
if hlc>(Dc-Dt)/2 € hlc<(Dc+Dt)/2 $€make sure hlc in correct 
range 
hw-hlc- (Dc-Dt)/2; 
else 
hw=-1;%tube is plugged or empty 
end 
% 
if hw~=-1 Stube not plugged 
phi w = acos(1-2*hw/Dt); 
alpha b = 1-(phi w - sin(phi w)*cos(phi w))/pi; 
else 
alpha b 
end 
$alpha b is the vapor volume fraction at the weir 


1; 


$use eqn () 

mdotw - 

rhol* (Dt/2)^2*pi*(1-alpha b)*sqrt (g*Dt*cos (theta) /a1); 
end 
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Bonction [phi c]=interp(alpha c) 

3 = 

Seumetion [phi c]-2interp(alpha c) 

N = 

function that finds phi given alpha using bisection search 


phi c-pi/2; 
bottom=0; 
top=p1; 
diff=1; 
while abs(diff)»5e-9, 
ENUL p5ha c-(l1—-(phi c-sin(phi c)*cos(phi c))/p1i); 
E رد‎ > 0 5 
bottom-phi c; 
pENEC-phi ct(top-phi c)/2; 
else 5 
top-phi c; 
Pepa c-(phi c-bottom) /2; 
end 
end 
end 
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function [j1,jv]= findjvjl(hl,Pprime,h, press, jvdat, jldat) 
% 


$This macro does a 2 way interpolation on jv and jl data 
compiled $as a function of both hl and P’ 

$hl is liquid height in the drain pipe 

$Pprime is dP/dx for the tube 

$press and givenh are the P’ and h data used to compute the 
jvdat %and jldat tables 

$ 

+function [jl,jv]2 findjvjl(hl,Pprime,h,press, jvdat, jldat) 
$ 


$find out where hl and Pprime fall in the range of h's and 
$press's available 


% 


$note that h and press contain all the possible values of hl 
and %Pprime 

% 

a=find(hl<=h & hl>=h(1)); %first element of a is the first h 
Shigher than given hl 

b=find(Pprime>=press & Pprime«-press(1)); 

% 


$now select the 4 jv/jl values that bracket the hl and P' 
given 

% 

if length(a)~=0 & length (b)-= 0 %hl and Pprime were in the 
$correct ranges 


% 
$fix pressure index (at lower bracket value, one less than 
5l) 
if ط‎ )1( << 
۱ / $in case pressure index already at lowest 
$value,don't decrement 
else 
b -b(1)-1;$fix pressure index at next lower value 
end 
$drop hl index to lower bracket 
if a(1)-2- 
ISE) Fl; $so a -2 will be 0,not negative 
else 
a =a(1); 
end 


jvhl-jvdat(41*(a -2)*b );$1ower bracket hl and Pprime data 
jlhl-jldat (41*(a (+ ) i 


$bump hl index back to upper bracket 


if b(1)-2- 
b -1; $in case pressure index already at lowest 
$value,don’t decrement 
else 


b -b(1)-1;$keep pressure index at next lower value 
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end 

EN a(1); 

jvh2-jvdat(41*(a -1)*b );$upper bracket hl, lower bracket 
$Pprime data 

۳۶۸ 11 03۲ )41* ج)‎ -1)-4b ); 


sbump pressure to upper bracket 


$fix pressure index (at upper bracket) 


if b(1)== 
f -2; Sin case pressure index already at lowest 
$value,bump up 
else 
f -b(1); %tput pressure index at upper bracket 
end 


$drop hl index to lower bracket 
if 2 )1( << 
g_=a(1)+1; 
else 
g_=a (1); 
end 
jvh32jvdat(41*(g -2)-f );$10wer bracket hl, upper bracket 
P' 
jlh3=jldat (41*(g -2)+f _); 


$bump hl index back to upper bracket 


if b(1)-2- 
f -2; $in case pressure index already at lowest 
$value,bump up 
else 
Zo), put pressure index at upper bracket 
end 
g_=a(1); , 


ve mica 41۰ )9 —1)+f£ );*upper bracket hl and P" 
jlh4=jldat (41* (g_ =1)+£_ yrs 
% jvh4-jvdat(11*(g -1)4f ); supper bracket hl and P’ 
$  jlh4- jidat(11*(g -1)-4f ); 


$interpolate between h values 
if 2 )1( << 
hl=h (1); 
else 
h1=h (a (1) -1); 
end 
h2=h (a(1)); 
h3=h1; 
házh2; 
Sthese if’s are required since hl might = h2 and the same 
for h3 52۳80 4 
if h2== 
jv5=jvhl; 
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else 

E-(brI-hl)./(h2-h1l)* (jvh2-jvh1)-4jvh1; 
end 
if h2==hl 

jv6=jvh3; 
else 

jv6=(h1l-h3)./ (h4-h3) * (jvh4-Jvh3) +jvh3; 
end 
if h2==hl 

:۳ 5 11 
else 

۲۱۳۰ ۱1۳1 -1:1( ۰7 )2-۳1( * )31۳27-11۳1( +311 
end 
if h2==h1 

711 62 7 3 
else 

Nos (hI-h3)./(h4-h3)* (j1h4-j1h3) *tj1h3; 
end 


<interpolate between 5 6 for Pprime 
1f b(1)== 
pl=press (b (1)); 
else 
pl=press (b (1) -1); 
end 
1f b(1)== 
p2=press (b (1)+1); 
else 
p2=press (b (1)); 
end 


Dr (Perime-p1)./ (p2-p1)* (Jv6-jv5)+3v5; 
M PEor:ime-pl)./(p2-p1)*(j16-j15)-47315; 


else 
b 
a 


disp('hl and Pprime were not in the correct ranges for... 


۳1 ۳ ۱ 1 ۰10۰ ( 
end 
end 
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